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Abstract. Strongly anisotropic critical systems are considered in a d- 
dimensional film geometry. Such systems involve two (or more) distinct 
correlation lengths £^ and £ a that scale as nontrivial powers of each other, i.e. 
£,a ~ So with anisotropy index S / 1. Thus two fundamental orientations, 
perpendicular (_L) and parallel (||), for which the surface normal is oriented along 
an a- and /3-direction, respectively, must be distinguished. The confinement of 
critical fluctuations caused by the film's boundary planes is shown to induce 
effective forces Tc that decay as JFq oc — (d/dL)Aj_ ^ L~ ^-'H as the film thickness 
L becomes large, where the proportionality constants involve nonuniversal 
amplitudes. The decay exponents Cj.,11 an d the Casimir amplitudes Ajii are 
universal but depend on the type of orientation. To corroborate these findings, 
n- vector models with an m-axial bulk Lifshitz point are investigated by means of 
RG methods below the upper critical dimension d*(m) = 4 + m/2 under periodic 
boundary conditions (PBC) and free BC of an asymptotic form pertaining to 
the respective ordinary surface transitions. The exponents II are determined, 
and explicit results to one- or two-loop order are presented for several Casimir 
amplitudes A^|. The large-n limits of the Casimir amplitudes A^ c /n for 
periodic and Dirichlet BC are shown to be proportional to their critical-point 
analogues at dimension d — m/1. The limiting values Aj^jpoo ~ lim n _ >oc II 
are determined exactly for the uniaxial cases (d, m) = (3, 1) under periodic BC. 
Unlike A^^, the amplitude A^ B ^ is positive, so that the corresponding Casimir 
force is repulsive. 
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1. Introduction 

When critical fluctuations in a medium are confined by the presence of macroscopic 
bodies such as walls, long-range effective forces between these bodies are induced. This 
phenomenon, first pointed out by Fisher and de Gennes [1], is the thermal analogue 
of the familiar Casimir effect between metallic objects caused by fluctuations of the 
electromagnetic field [SJ. It has attracted considerable attention during the past 15 
years [3J 13] • Its indirect experimental verification through the thinning of 4 He wetting 
layers near the lambda transition [SJ O [7] some time ago and the more recent first 
direct measurement of such Casimir forces for binary fluid mixtures [8, 9,110] are likely 
to spur further interest and increasing research activities in this field. 
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Previous studies of critical Casimir forces in statistical physics have focused 
exclusively on macroscopic media exhibiting isotropic scale invariance in the absence 
of confining walls and macroscopic bodies. The characteristic feature of such systems 
is that they become self-similar when distances Ax along arbitrary directions are 
rescaled by a scale factor I = Ax /Ax'. Alternatively, one can say that the bulk 
correlation lengths characterizing the decay of correlations along all directions exhibit 
the same power-law divergence ~ \T/T C — as the temperature T approaches the 
bulk critical temperature T c . 

In this paper we shall be concerned with strongly anisotropic scale-invariant 
systems. In bulk systems of this kind there exist one or several principal directions 
along which coordinate separations Ax a must be rescaled by a nontrivial power i 6 of 
the scale factor i = Axp/ ' Ax'p associated with the remaining principal directions in 
order to have self-similarity of the initial and transformed systems. The asymptotic 
power laws 

Up-\t/t c -\\-^ (l.i) 

of the respective bulk correlation lengths £ Q and £p are governed by distinct exponents 
v a = Ovp and i/p. Such systems are ubiquitous in nature. Important examples of 
equilibrium systems whose (multi) critical equilibrium states exhibit anisotropic scale 
invariance ( ASI) are systems at Lifshitz points [TTJ Q21 [13] and liquid crystals [14] [15] . 
Furthermore, ASI is a common feature of the stationary states of many nonequilibrium 
systems [15]. We shall show that the Casimir forces induced by confinement of 
fluctuations in systems exhibiting ASI differ qualitatively and quantitatively from their 
analogues for isotropic scale-invariant systems. 

Consider a system confined by two parallel planes at a distance L, each of which 
has area A. The reduced free energy of the system per area A can be written as 

» L/ b (T, . . .) + / S (T, . . .) + / res (L; T, . . .) (1.2) 

in the limit A — > oo at fixed L, where /b and / s are L-independent bulk and surface 
excess densities. The ellipsis in fh(T, . . .) stands for additional thermodynamic bulk 
fields (such as magnetic field or chemical potential); the one in f s (T, . . .) represents 
both bulk fields of this kind as well as additional surface variables (such as surface 
interaction constants). The L-dependence resides in the so-defined residual free-energy 
density / res (L;T, . . .). For given medium and boundary planes, this function decays 
at the bulk critical point as 

fr eS ,0Tit( L ) « A BC A^ (i/Ll)" C (1-3) 

as L — > oo, where A\ and L\ are metric factors (units of area and length, respectively). 
At conventional critical points exhibiting isotropic scale invariance, one has £ = d — 1 
for the decay exponent. Furthermore, the metric factors are chosen to satisfy 
A\ = L^ 1 for systems whose correlation regimes are hyperspherical (i.e., whose bulk 
correlation lengths £ w £q\T/T c — l\~ v characterizing the decay of correlations along 
different spatial directions diverge with the same exponent and have equal amplitudes 
£o)- This choice guarantees that A\ and L\ drop out of equation (|1.3|) . 

The "Casimir amplitude" A BC depends on gross properties of the medium 
(universality class) and the boundaries, namely, the boundary conditions (BC) that 
hold on large length scales [T7] , but is independent of microscopic details (universal) . 
To avoid confusion, let us explain how the metric factors A\ and L\ ought to be 
chosen when the hypersphericity condition on the correlation regime is violated in 
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the way weakly anisotropic critical systems do [TS]- For those, the bulk correlation 
lengths £ associated with distinct directions diverge with the same critical exponent 
v but involve several different amplitudes £o, so that their correlation regime is 
hyperellipsoidal. Their critical behavior can be expressed in terms of that of 
standard isotropic systems such as the conventional 4> 4 model [HI [19] . The required 
transformation that makes the correlation regime hyperspherical — a similarity 
transformation followed by a rescaling of the principal radii — generally changes A\ 
and L\ into transformed values A\ and L[ (see p. 15-17 of [19]]). To define A BC 
for weakly anisotropic systems in a consistent manner, one must simply choose 
and L\ according to our rule for isotropic systems, requiring them to drop out of the 
corresponding analogue of equation (11.31) . 

The anisotropy one encounters in weakly anisotropic critical systems is of a 
fairly harmless kind: it can be transformed away, absorbed by a proper choice 
of (nonuniversal) coordinates. This is not the case for systems exhibiting ASI. 
Their anisotropy is of a more profound type. This has important consequences for 
fluctuation-induced forces. General aspects of the orientation of the boundary planes 
matter. Two fundamentally distinct orientations must be distinguished: parallel (||), 
for which the normals to the boundary planes are oriented along a /3-direction, and 
perpendicular (_L), for which the normals are parallel to an a- but perpendicular to 
all /3-directions. If one has m a- and d — m /3-directions in a (i-dimensional system, 
then the corresponding decay exponents are given by 

C|| = d - m + 9 m - 1 and £j_ = (d — m)/9 + m — 1 . (1.4) 

In the cases of m-axial Lifshitz points, which we explicitly consider in the following, 
the value of the anisotropy exponent 9 is close to 1/2; one has 9 = l/2 + 0(e 2 ), [20l[2Tj 
where e = d* (m) — d is the deviation of the bulk dimension d from the corresponding 
upper critical dimension d*(m) = 4 + m/2 of the system. 

The associated Casimir amplitudes A^ c and A 1 ] 50 depend on the BC, and 
generally differ for perpendicular and parallel orientations. How they should be 
defined in order to avoid trivial dependences on nonuniversal metric factors needs 
to be clarified. The results f| 1 .4|) for the decay exponents can be understood by simple 
scaling arguments. The residual reduced free energy density f Tes (L) has dimension 
1/A where A scales as C™^ - ™ -1 and £™ _1 £^ _m for parallel and perpendicular surface 
orientations, respectively. If £ a and £^ are the only relevant lengths besides L, then 
equation (11.41) should holdJH 

To substantiate these claims and verify explicitly that the Casimir amplitudes 
differ for parallel and perpendicular slab orientations, we shall investigate a familiar 
class of prototype n- vector models exhibiting ASI — namely, 0(n) <p A models with an 
m-axial bulk Lifshitz point (LP) [TTJ [T^J [T3] in a slab geometry. The cases of parallel 
and perpendicular orientations of the boundary planes will both be studied under 
periodic (PBC) and free (FBC) boundary conditions. However, when considering 
FBC, we shall restrict ourselves in two ways: We assume that the BC that result 
in the large length-scale limit (i) do not break the 0{n) symmetry and (ii) arc 
associated with the respective most stable renormalization-group (RG) fixed point. 
For parallel orientation this simply means that Dirichlet BC <p = hold asymptotically 
[22l [23l [24l [25] . In the case of perpendicular orientation, two BC hold on either 

X Note that since £ a and £^ are bulk correlation lengths, they diverge at the bulk critical point. In 
later sections we shall also consider finite-size correlation lengths. These remain finite at the bulk 
critical point when L < oo. 
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boundary plane. These simplify in the large-length-scale limit to the conditions that 
both the order parameter <fi and its normal derivative d n <j> vanish. 

In the next section we introduce the models and specify their Hamiltonians 
including the boundary terms they involve in the cases of free boundary planes. In 
section [3] we first give the mesoscopic BC that result from the boundary contributions 
to the classical equations of motion in the cases of parallel and perpendicular 
slab orientations. Assuming that the values of the surface interaction constants 
comply with the above-mentioned conditions (i) and (ii), we then recapitulate which 
asymptotic boundary conditions are encountered in the limit of large length scales. 
In section [4] we recall the background on the renormalization of these models at and 
below their upper critical dimensions d*(m) = 4 + m/2 and the field-theoretic RG 
approach to bulk and surface critical behavior at LP required for our subsequent 
analysis of the Casimir forces. We then turn to the calculation of fluctuation-induced 
forces at the LP. The case of parallel slab orientation is dealt with in section 15.11 that 
of perpendicular orientation in section 15.21 Section |B] provides a brief summary and 
concluding remarks. Finally, there are 3 appendixes describing technical details. 

2. Models 

The models we consider involve an n-component order-parameter field 4>{x) = 
(<j) a (x),a = 1, . . . ,n) defined on the slab 03 = R d_1 x [0,L] of e?-dimensional space 
R d . We write position vectors as x = (y,z), where y € R d_1 and z £ [0,L] are 
the coordinates alongside and across the slab, respectively (see figure [TJ. Parallel 
orientation means that z is a /3-coordinate, perpendicular that it is an a-coordinate. 
Without loss of generality, we can take the first m Euclidean axes as a-directions. We 
choose z to be the j3- or a-coordinate with the largest index, so that z — x<i and z = x m 
for parallel and perpendicular orientations of the boundary planes, respectively. 

The slab is assumed to have infinite area A = oo of its boundary planes. We 
can therefore choose PBC along the corresponding d — 1 principal y-directions for 
convenience. Depending on whether we take PBC or FBC in the z-direction, the slab 
03 has no boundary, c>03 = 0, or its boundary consists of the two (d — l)-dimensional 
confining hyperplanes 03 1 at z = and 032 at z = L. In the latter case, we orient the 
boundary such that the normal n on 903 = 03 = 03i U 032 points into the interior of 
03. 




Figure 1. Slab geometries considered: For perpendicular and parallel 
orientations of the surface planes at z = and z = L, the z-axis is along the 
o-direction x m and the /3-direction z^, respectively. The remaining y-dircctions 
are a- or /3-directions. Further, n is the inward-pointing normal. 

Ignoring any long-range interactions, we choose the Hamiltonians to be of the 
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form 



2 

U = \ £ b ( : r)dl/ + ( 5BC,frccV / £ 



(x) dA , 



(2.1) 



where /^(^O an d Cj(x) are functions of the field 4>{x) and its spatial first and second 
derivatives d 1 (f>(x) = dcj>{x)/dx 1 and d 1 d 1 >4){x). Following [20] ET], we choose the 
bulk density appropriate for the description of the (multi) critical behavior at m-axial 
Lifshitz points as 

d 



4! 



(2.2) 



a—l /3=m+l 

On the level of Landau theory, the LP is located at f = p = 0. 

When m > 1, this choice of Ch ignores potential anisotropies breaking the 
rotational invariance in the m-dimensional Euclidean a-space R m , such as a term 
proportional to the hypercubic invariant J2™=i(da4>) 2 ■ According to the two-loop 
calculation reported in [27], such a contribution is relevant in the RG sense. We omit 
it here, as well as similar less symmetric terms, for the sake of simplicity. In the 
uniaxial case m = 1, this is no restriction. 

Throughout the paper, we shall consider PBC and FBC. In the former case, 

0(y,z + L) = 0(y,z) (2.3) 

and the boundary terms in equation (|2.1I) are absent, as indicated. 

In the case of FBC, different boundary densities C^(x) and Cf{x) (dictated by 
relevance/irrelevance considerations) must be chosen to define appropriate minimal 
models for slabs with parallel or perpendicular orientation. Work on semi-infinite 
system^ [22 O EH EH [26] suggests the choices 



and 



C ]{x) = -c 3 4? + 



rn— 1 

+ £ 

a=l 



(2.4) 



b J <f>d n <t> + -\f(d n <j>Y 



■ A ! (d a (j>) 2 + fj {d a 4>)d n d a 4> 



(2.5) 



where we temporarily allow all bare interaction constants &j, . . . , fj to take different 
values on the two boundary planes QSj, j = 1,2. 



3. Mesoscopic and asymptotic boundary conditions 

As discussed elsewhere [28[ [29], for actions of the form (|2.ip . the boundary 
contributions to the classical equations of motion give us "mesoscopic BC" that hold 
in an operator sense. We call them mesoscopic because they hold on the length scale 
beyond which the chosen continuum description is valid. They must be distinguished 
from the asymptotic large-length-scale BC one encounters at criticality. Since the 

§ For general background on the field-theoretic RG analysis of critical behavior in semi-infinite 

systems, sec 28 and |29l . 
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boundary contributions to the action differ in the cases of parallel and perpendicular 
orientations, so do the mesoscopic BC. We first consider the case of parallel slab 
orientation. 



3.1. Parallel slab orientation 

In this case, the action is defined by equations (|2.1[) . (12.21) and (|2.4[) . The mesoscopic 
BC become [H [23] 

(3.1) 



d n 4>{x) = U - Xj jr gA 0( X ) , xexj 

^ Q = l ' 



The bare interaction constants Cj and have engineering dimension [cj] = fi 
and [Xj] = where /i is a momentum scale. The former are relevant at Gaussian 
fixed points with Cj = Xj =0. Beyond Landau theory, the associated renormalized 
quantities Cj, which are proportional to the bare deviations 5c j = Cj — c sp from a 
cutoff-dependent special value c sp , take over the roles of relevant surface scaling fields 



[24 UJ We assume that the initial values of the renormalized surface variables Cj and 
Xj lie in the basin of attraction of the fixed point with Cj = oo describing the so-called 
ordinary surface transitions of the semi-infinite systems with one surface plane at *Bj 
and the respective other at z = oo or z = — oo. As is known from |22) . |23j and 
[2~I] . Dirichlet BC hold at this fixed point at both boundary planes 03 j. To study the 
corresponding asymptotic behavior, one can set the bare values Cj = oo. Then the 
regularized bare theory satisfies Dirichlet BC. The values of the bare variables Xj > 
do not matter. To investigate the critical Casimir forces in this case of parallel slab 
orientation and large-scale Dirichlet BC, we can simply set b\ — 62 — 00 and drop the 
boundary contributions oc Xj. Equivalently, we can impose Dirichlet BC 

0(x)=O, VzeSiUQSi, (3.2) 

on the regularized bare theory and drop the boundary terms Cj dA of the action. 

3.2. Perpendicular slab orientation 

In this case, the action is defined by equations (|2.ip . (|2.2[) , and (|2.5p . Since z now 
is an a-direction, the classical equation of motion is of fourth rather than of second 
order in d z . In conformity with this we find two (instead of one) mesoscopic BC at 
either boundary plane, namely |25l 126] 

^ad 3 n + (bj - p)d n + cf - [X [ \ + (fj - a)d n ] £ dl^ix) = (3.3) 

and 

m—l 

- & dl + Xfd n + bj - (f) +cr)^2dl 4>(x) = for x e . (3.4) 

a=l 

The problem of boundary critical behavior at Lifshitz points has been investigated 
to a much lesser degree for the case of perpendicular orientation of boundary planes. 

I When deviations p oc p — plp of the bare bulk variable p from its value plp at the bulk LP are 
considered, the corresponding linear surface scaling fields g Cj actually are linear combinations of Cj 
and p, the renormalized counterpart of p. Details, which will not be needed below, can be found in 

EH. 
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However, as shown in 25 and [26 , one can benefit from analogous simplifications 
when analyzing the asymptotic behavior at the corresponding so-called "ordinary" 
surface transition (described by the most stable fixed point with bulk LP behavior). 
Let us recapitulate the essence of the argument leading to the conclusion that 

<p(x) = 0, d n (j>{x) = 0, \fx e 05 1 U 03 2 , (3.5) 

are the appropriate large-length-scale BC to consider. Noting that d n now has 
engineering dimension /i 1 / 2 , one finds that the interaction constants of the boundary 
densities Cj~ with the largest momentum dimensions are cj~ ~ /i 3 / 2 and bj ~ /i. Thus 
both give rise to scaling fields that are relevant at the Gaussian fixed points where 
they vanish. For generic positive initial values, they will be driven to oo under RG 
transformations p — > fj£ in the infrared (IR) limit £ 0. We therefore take the limits 
Cj —> oo and bj — > oo at the outset. Upon dividing the boundary conditions (|3.3p 

and p.4[) by cj- and bj , respectively, and scaling all interaction constants according to 
their p dimensions, one arrives at the asymptotic conditions (I3.5P for the regularized 
bare theory. 

The upshot of these considerations and those in the previous subsection is the 
following. To investigate the asymptotic form of the Casimir forces at the LP in 
the case where the surface interaction constants are subcritically enhanced at both 
boundary planes, we can impose the Dirichlet BC (13.2[) and the BC (13.5[) . depending 
on whether the slab orientation is parallel or perpendicular, and omit the boundary 
terms J2j I<s *S ^ of the respective actions (12.11) . 



4. Background and renormalization group 

We wish to investigate the residual free energy / rcs in d = d* (to) — e dimensions by 
means of the RG approach. Our aims are to confirm the asymptotic behavior ()1.3[) and 
compute the Casimir amplitudes A^ f° r parallel and perpendicular slab orientation 
under the mentioned BC using RG-improved perturbation theory. To this end we 
must first supply some background on free propagators and the renormalization of the 
models. We begin by recalling a number of required results of the free bulk theory. 



4--1- Free bulk propagators 

For notational convenience we decompose the position vector x = (y, z) into its Tri- 
dimensional a-component r = (x a ) and (d — TO,)-dimensional /3-component s = (xp), 
writing x — (r, s). Here and below it is understood that a runs from 1 to m, while 
/3 = m + 1, . . . , d. Likewise, we use the notation q = (k,p) with k = (q a ) and p = (qp) 
for d-dimensional wave vectors q conjugate to x and their a- and /3-components. 

The free bulk propagator G^'"^ of the disordered phase in d dimensions follows 
from the Gaussian part of the bulk density £b given in equation (|2.2I) . For general 
non-negative values of p and f , it can be written as 

{d) p i(rk+sp 



where 



G^ m >(x\a;r,p)= , (4 .1) 



( d ) r Ad„ Am) Ad— m) 



(27T) 



(4.2) 
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denotes a conveniently normalized momentum integral. By rotational invariance in 
the r- and s-subspaces, the position dependence of G^' m \x\a; r, p) reduces to a 
dependence on the lengths r = \r\ and s = \s\ of the a- and /3-components of x. 

We shall use dimensional regularization. Unless stated otherwise, all momentum 
integrations will therefore not be restricted by a large-momentum cutoff A. Being 
interested in properties at the LP (such as the Casimir amplitudes A BC ), we will 
normally set p = r = and work with the free LP propagator G^f (x) = 
G^' m) (x\&;0,0). This can be written as [20JET] 

G { h d ' m \x)=a- m / i s- 2+ ^ m Av), v^a-^rs- 1 / 2 (4.3) 
with the scaling function 

**,<(«) = / fl -jfrntr, (4-4) 

where s is an arbitrarily oriented unit vector in R d ~ m . For general values of to and d, 
the latter function is a difference of generalized hypergeometric functions 1F2 (equal 
to a Fox- Wright "J function [21]). Its explicit form can be found, for example, in 
equation (A2) of [20] or equations (11) and (13) of [21]. For our purposes here, it will 
be sufficient to know its Taylor expansion 

mA ) ^ 2 2i Z! r[(2 + m + 2Z)/2) [ ' 

and its asymptotic expansion for large v, 



(-1)' 2 4 'r(2-6 + 2Q 
^ II r[(m-2 + 2e-4Z)/4] 



« 2 J - m -^7r^- - m ^^V^ ^ ^ , ^77T ^~ 4 '- (4.6) 

From the former one can infer the value at v = 0: 

$ m>d(0 ) = y 1 ~ £/2) . (4.7) 

T[(2 + m)/4J 

We shall also need the leading asymptotic behavior of & m ,d for large v below. 
According to equation (14.61) . we have 

$ m>d (v) ^ ^ £ " 4 + 0(v 2 '~ s ) (4.8) 

with 

.(oo) _ 2 3-m-2e (2e-6-m)/4 F ( 2 ~ £ ) f 4 

r[(m-2 + 2e)/4] ' 1 J 

Note that the coefficient vanishes on the whole line <i = 3 + to. This reflects the 
fact that the function 4> m ,3+m(u) decays ~ e~ v2 / A (cf. equation (19) of [2"T]). 



^.1?. Reparametrizations and renormalization group 

Let us introduce the (dimensionally regularized) Appoint cumulants (connected 
correlation functions) 

.Jxx, ...,**) = (0 ai • • • </w(^)) cum (4.10) 

and denote their bulk (infinite-space) analogues as G^J a N i x i' ■ • ■ 1 x n)- The 
renormalization of the latter functions in d — d* (to) — e < d* (to) dimensions has 
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been explained in detail in 123 H2| ■ Their ultraviolet (UV) divergences induced 
by the singularity of the free bulk propagator G^' m \x) at x = can be absorbed by 
means of the reparametrizations 

<p = Z]J 2 <pn, a = Z a a , u a~ m/A F m ^ = pf Z u u , 

T-TLP = fl 2 Z T ( T + A T p 2 ) , (p- p L p) CT' 1 ' 2 = HZ pP , (4.11) 

where p is a momentum scale. Following JT] and [22], we choose the factor F mj£ that 
is absorbed in the renormalized coupling constant as 

F r(i + e /2)r 2 (i- £ /2)r(m/4) 

m,e (4^)(8+'»- 2£ )/ 4 r(2-e)r(m/2) ' l ' J 

The LP is located at (f, p) — (tlp,/5lp)- In a theory regularized by means of a 
large-momentum cutoff A, the renormalization functions tlp and /?lp would diverge 
~ A 2 and ~ A, respectively. In our perturbative approach based on dimensional 
regularization, they vanish. Results to order u 2 for the renormalization factors 
Z . Z p , Z T and Z u may be found in equations (40)-(50) of [3T]. The function A T is 
given to 0(u) in equation (17) of [2"2"] . 

Upon varying /i at fixed values of the bare variables it, a, p, and f , we find that 
the renormalized bulk functions G]^r = N ^ 2 Q^' satisfy the RG equations 



N 



V.+ '-m) a„ ;R (*l> ■ • • • = ' - U •» 1 

with 

©„ = A*^+ E Mr ( 4 - 14 ) 

The beta functions j3 g are defined by 

P g = pd^og , g = u,a,p,T, (4.15) 

where c^jo means a derivative at fixed bare interactions constants. They can be 
expressed in terms of the exponent functions 

r) g (u) = pdplolnZg , g = <j>,u,<T,p,T, (4.16) 



and 



as 



b T (u) = A T [p df,\ In A T +n T - 2n p ] (4.17) 



-u[e + T] u (u)}, g 



u. 



-p[l + r]p{u)], g = p, 



-t[2 + r) T (u)] - b T (u)p 2 , g = T. 

The reparametrizations (|4.11[) also suffice to absorb the UV singularities of the 
iV-point functions Gj^2 a for films of finite thickness L under PBC, irrespective 
of whether the orientation is parallel or perpendicular. This can be concluded from 
the form of the free film propagator (see, e.g., [28j chapter IV]), 



GS°l(*-*0= E G^ m \x-x'+jLz), (4.19) 
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where z is a unit vector along the z-axis. The crux of the argument is that all primitive 
UV singularities must be induced by the j = summand (the bulk term) , because the 
remaining j 7^ summands are finite at x = 0. Accordingly, the RG equations (|4.13|) 
hold also for the film cumulants Gi N } „ under PBC. 

In the case of FBC one generally expects additional primitive UV singularities 
localized at the boundaries. Whenever they occur, counterterms with support on 
the boundary planes 53^ ("surface counterterms") must be added to the action. 
However, we can benefit from the simplifying features of the chosen Dirichlet 
BC (|3.2p and (13.51) for parallel and perpendicular orientations. These BC ensure 
that the renormalized correlation functions G\ N ) „ = Z, N ^ 2 g{ N ) , defined via the 
reparametrizations (|4.11[) become UV finite even when some (or all) fields 4> ai (xi) 
are located on the boundary. Surface counterterms would be required in the case 
of parallel slab orientation if we wanted to renormalize cumulants involving normal 
derivatives d n 4> a in addition to fields 4> a {x) with x e [221 For perpendicular 
surface orientation, such cumulants would still be UV finite owing to the stronger BC 
that both <fi and d n <p vanish. On the other hand, cumulants involving the boundary 
operators d^cj) would require surface counterterms |26) . 

In order to investigate the residual free energy f Tes (L; . . .) by means of the RG 
approach, we must know how to renormalize this quantity. The bulk and surface 
excess free energy densities /b and / s are known to require additional additive 
counterterms. The subtractions these counterterms provide do not depend on L. 
Thus, they cancel in the difference of the free-energy densities defining f res (L; . . .). 
Using the reparametrizations (|4.11j) to express bare variables in terms of renormalized 
ones therefore gives us UV finite renormalized residual free-energy densities 

/rcs;R.(£; T, p, U, (7, fl) = f ICS (L; T , f), U, &), (4.20) 

both for PBC as well as either one of the Dirichlet BC (pT2")) and (|53)l . 
As an immediate consequence we obtain the RG equation 

T> /i/res;H.(£; T, p, U, (7, ft) = . (4.21) 

Solving it at the LP r = p = via characteristics and using dimensional analysis gives 
/ rcs;R (L; 0, 0, u, a, fi) = / rcs;R [i; 0, 0, u(£), a{£), fit) 

,2ff2\ m / 4 



G^)*-™- 1 ( ) U-A^ld; 0, 0, u(£), 1, 1] 



\d-r 



2i>2 \ (m-l)/4 



(/"■i" m [j^) / ros;R [L(^ 2 M^)) 1/4 ;0,o,^),i,i], s 



(4.22) 



where the variable <^ specifies the orientation. Further, u(£) and o~(£) are examples of 
running variables g(£) defined in the usual way as solutions to the flow equations 

e^m = P B [fV),p(£)Mt)MQ], a = *,..., a, (4.23) 

subject to the initial conditions 

g(£ = 1) = g. (4.24) 

In the large- length-scale limit £ — > 0, the running coupling constant approaches 
the IR-stable root u* , whose e expansion 

u* = -^-+0(6 2 ) (4.25) 
n + 8 
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is known to second order |21) but will be needed only to first order in the following. 
The variable a(£) behaves as 

ail) ps E*(u)i- r '^ u *K, (4.26) 

where the fixed-point value of r\ a is related to the anisotropy exponent 9 by 

6= [2 + 7 ?CT ( M *)]/4 (4.27) 
and E*(u) = E a (u* ,u) means the value of the trajectory integral 

^ (u ' w) = exp U du w I 

= exp ^ ^ {40 - 2 - vMm) (4-28) 

at u = u* . 

We now choose I such that the scaled L-dependent arguments of / res ;R m 
equation (|4.22p become 1, obtaining 

f ^d" 20 )/ 2 (£»-™/ 4 Af c L^ii , ? = ||, 

L^Too 1 (d-m)(2fl-l)/(2fl) ^* (T )(d-m)/(40) A^ C L _Ci , C = _L. 

(4.29) 



/SjR(i;0,o,«,o-,A*) 



with 

Af C = /l C ; R (l;0,0,u*,l,l), (4.30) 

where are the decay exponents introduced in equation (|1.4j) . BC stands for either 
PBC or else the asymptotic ones associated with the ordinary surface transitions, 
at both boundary planes. Depending on the orientation, the latter are given by the 
Dirichlet BC (I3.2[) (c = ||) or equation Q3.5P (<; =-L). In either case, we use the notation 
BC = (O, O) to refer to these boundary conditions. 

4-3. Casimir amplitudes as universal ratios 

Note that the result (I4.29[) involves, besides the momentum scale [i and the variable a 
(both of which are there for dimensional reasons), the nonuniversal amplitude E*(u). 
On the other hand, it is not difficult to see that the amplitudes A^ c can be written 
as universal ratios of nonuniversal quantities. To this end, let us introduce the second 
moment finite-size correlation lengths 

-. m 

^ = -E(^-^) 2 ( 4 - 31 ) 

a=l 

and 

&**=4Z^ E (** " 4) 2 ' ( 4 - 32 ) 

/3=rn+l 



where the averages h(x,x') are defined by 



[„ d d x L, d<V h(x, x') G% (as, x') 

h(x, x>) = Jv - J * - ^ ^ £3 \; . (4.33) 



I<a ^ dx Jsn d dx ' @l--r( x > x> ) 
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As indicated, these quantities depend on the orientation <; — | or <; =_L whenever 
L < oo. In the limit L — > oo they approach the corresponding (orientation- 
independent) second-moment bulk correlation lengths £b;a = £oo ; a,? and £b ;/ 3 = 
^oo ; /3.<;- Note that we assumed the bulk density (|2.2[) to be isotropic in both the 
a- and /3-subspaces. Therefore, the bulk correlation lengths can equivalently be 
written as bulk moments £, 2 . a = [aJf] b and ^ ;/3 = [x 2 ] h for arbitrary choices 

of a = 1, ...,m and /3 = m+ 1, ...,d, where [/i(a;)J b denotes the bulk average 

J Rd d d xh(x) G^ K {x)/ J Rd d d x G h 2 ^(x). A violation of the bulk density's isotropy in 
/3-space would correspond to what is usually called "weak anisotropy" . It means that 
the bulk moments , associated with different values of P have the same power- 
law singularity ~ f~ Vla as ^ ne Lp j s approached along a path with p = but may 
have different (nonuniversal) amplitudes. Getting rid of such a weak anisotropy can 
be achieved in a straightforward fashion by transforming to order parameter densities 
that yield a squared gradient term in /3-space of the form assumed in equation (I2.2[) . 
To achieve this, one can proceed along lines analogous to those taken in (T8l [19] to 
cope with weak anisotropy at critical points. 

Anisotropies in a-space are more dangerous. The two-loop RG analysis [27] of 
isotropy-violating linear combinations of (d 2 <p / dx 2 a ) 2 to the bulk density Cb revealed 
that the fixed point in d*{m) — e dimensions is unstable to such perturbations. For 
this reason we explicitly rule out the presence of such isotropy-breaking terms. 

Now, the RG equations D ^ £l ;q , i? = T>^ £l-b <; = in conjunction with dimensional 
analysis yield 

£i;a, S (u, T, p, a, p) = (p£y 1/2 CT 1/4 il-a^U, f, p, 1, 1) 

u-V [E:(u)afU- 6 ^ a Ju*,f,p,lA) (4.34) 



e-to 



and 



where 



W%,c(«'>np>l.l), (4-35) 



f pih for? = ||, 

L= [ M Va^/a^-i/4 £ a y^H e L for t =i_. (4 ' 36) 

Setting t = p = and choosing the flow parameter 1 = 1^ again such that L = 1 
gives 

f = A 172 (£*^) 1/4 ^, i e for c = ||, 

< !C ; (M,0,^) r « <| „^ _^ BC ; ■ „ " (4.37) 



L = ££l L for C =i_, 



and 

^(u.O.O.a,/*) « rlH " , . rl „ ,. ." <4.oM 



H| j ( jji=e^ii fo^ 



L-oo ^ S BG L l/fl = ^^^-^ gBC ± L Vfl for ? =±j 

where we have introduced the universal numbers C*™^ = £1^7 c( u *! 0, 0, 1, 1) with 

7 = a, j3 and c = || , I . Owing to the presumed isotropy of our model in both the 

a- and /3-subspaces, the amplitudes 3*^ are the same for all possible m choices of a. 
Likewise, 3^ is independent of the choice of /3. 
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Using the above results and definitions, one easily concludes that the Casimir 
amplitudes can be written as 

Af c = lim f™ R (L;0,0,u,o-,ri 



ggk>,0,0,^) 



0,0,<7, At) 



-i d—m—Sr 



*BC 



(4.39) 



i.e., the dependences on /i, er and the nonuniversal scale factor drop out of these 
ratios. Note that this definition would even work if weak "diagonal" anisotropics were 
allowed in both the a- and /3-subspace. By such anisotropics we mean the kind that 
can be transformed away by simple rescalings of the Cartesian coordinates x a and xp, 
respectively. In other words, we assume that the metrics in the a- and /3-subspaces 
provided by the respective gradient square terms in the Hamiltonian remain diagonal 
[15] . The case of lattice systems that lack hypercubic symmetries in the a- and /3- 
subspaces, such as systems with monoclinic lattices, requires somewhat more thought 
since the metric in full space must be diagonalized. To this end, the considerations 
made in [T5] Q15] for weakly anisotropic 4> 4 models near critical points (CP) must be 
appropriately adapted and generalized. 



5. Calculation of fluctuation-induced interactions at the bulk LP 

Having established the RG predictions for the asymptotic behavior of the residual 
free energy density, we now turn to the calculation of the Casimir amplitudes Af c 
by means of RG-improved perturbation theory in d = d* (m) — e dimensions. We first 
consider the case of parallel slab orientation. It is computationally less involved than 
that of perpendicular film orientation. Furthermore, one can proceed in close analogy 
to the previous analysis of Casimir forces at critical points [30] [23 GUI [23] ■ 



5.1. Parallel slab orientation 

We begin by considering the exactly solvable Gaussian case with u = 0. In this 
noninteracting case, we have a = a, f = /i 2 r, and p — pp. To obtain the Gaussian 
analogue of the result for / rcs , we must insert the classical value 1/2 for the anisotropy 
exponent 9 both in equation (|4.29l) and expression (ll.4[) for the decay exponent 0| • In 
addition, the nonuniversal amplitude E* must be dropped. This yields 

/ r B es c (i;0,0,0,a,At) = <7- m / 4 Af3 ( 5>1) 

In | Appendix A| we show that the Casimir amplitudes are given by 

Af G c (d,m,n) = 2 d -^A^\d, m ,n) 

r r[(rf-m/2)/2R(d-m/2) 

~ n0m _(d-m/2)/2 



with 



n {2-m)/4 



Cm ~ 2™r[(m + 2)/4]- (5 ' 3) 

The factor C m is a ratio of geometric factors one encounters in the calculation 
of Feynman integrals for bulk critical behavior at LP and critical points (CP) that 
involve the respective free bulk propagators G\f' m \o\a; 0, 0) and G^' ^ (0|<j; 0, 0) at 
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coinciding positions x' — x (tadpole graphs). Two examples of its occurrence are the 
ratios 

n _ F m . e _ $m,4+m/2-e(0) 

<^m — — T 7777 , 

F , e $0,4-e(0) 

where the &m,d(0) are the values of the functions (|4.4[) at u = 0. 

To understand the origin of C m , note that Gj^' m) (0|<7;0,0) according to 
equation (|4. 1[) is given by a momentum integral of the form J J^™" 1 /i(p 2 + 
crfc 4 ) with ft.(g 2 ) = 1/g 2 - Such integrals can be expressed in terms of the 
(d — m/2)-dimensional momentum integrals J^ d ™^ 2 ' > h(q 2 ): using hyperspherical 

coordinates in k and p space, one can transform to the variable K = er 1 / 2 fc 2 
and then introduce the polar coordinates p = qcosO and K = qs'md, following 
[34] . The integration measure p d ^ m ^ 1 fc™ -1 dpdfc of the radial integrations becomes 
£-m/4 (qcos0) d_ ' Tl-1 (gsin0)(" l ~ 2 ^ 2 gdgd#/2, and a straightforward calculation leads 
to 

(rn) p(d—m) p(d—m/2) 

/ h(p 2 +ak 4 )=a- m ^C m h(q 2 ). (5.5) 

k J p J q 

Returning to the above results (|5.ip and (|5.2[) . let us note that they must reduce 
to their known CP analogues in the limit m — > 0. The Casimir amplitudes A^p' q (d, n) 
of the free field theory with action density Cb = '}2™ = i{d 1 4>) 2 /2 can be read off from 
[351 equation (4.4)] and are explicitly given in [311 equation (5.6)]. Their analogues 
A^p^j(d, n) can be found in [2U equation (5.7), 3rd line]0 Taking into account that 
Co = 1, one sees that the corresponding expressions are indeed recovered from our 
equations (|5.ip and (|5.2p in the limit m — > 0. 

A straightforward consequence is that the Casimir amplitudes (|5.2[) can be 
expressed in terms of their CP analogues as 

Af£(d, m, n) = C m Ag£ G (d - m/2, n). (5.6) 

In fact, the analogue of this relation between A| c (d, m, n) and A|=5p(<i— m/2,n) turns 
out to hold in the limit n — ^ oo for 2 + m/2 < d < 4 + m/2. To make this statement 
precise, let us define the large-n amplitudes 

AJ?° (d,m) = Hm Af, (d, m, n)/n for BC = PBC, (O, O), (5.7) 

and their CP analogues A^p ^(D), where we assume that 2 < D = d — m/2 < 4, i.e. 
that D lies between the lower and upper critical bulk dimensions d*(m = 0) = 2 and 
d* (m = 0) = 4, respectively. Our claim is that these amplitudes obey the relation 

A|°(d,m) = C m Agg )00 (d-ro/2), BC = PBC, (O, O). (5.8) 

Its derivation is given in |Appcndix B| For PBC, this is straightforward. The 
case of (O, O) BC requires more thought because the loss of translation invariance 
perpendicular to the boundary planes entails that the exact solution in the limit 
n oo involves an effective z-dependent pair interaction, which must be determined 
self-consistently [3lfl 1371 13"5 ] . 

One may wonder whether analogues of the relations (|5.6p and (I5.8P might hold 
in general. We emphasize that there is no reason to expect this. The origin of their 
validity in the free-field case and the large-n limit is that the corresponding amplitudes 

% Note that a minus sign is missing in this equation of 1311 . This is an evident misprint. 
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involve only Feynman integrals of a sort that equation (|5.5[) can be applied0 However, 
at higher orders of the loop expansion Feynman graphs involving powers of free 
propagators between different points in position space occur. Relation (I5.5[) is not 
applicable to them. As we shall see below, the relations (|5.6[) and (I5.8[) carry over 
to the small-e expansion to the order of our calculation. The reason is the same 
as before: only Feynman graphs involving free propagators at coinciding points 
contribute. At order e 2 , three-loop graphs involving powers of free propagators 
between distinct positions contribute. These must be expected to invalidate the 
analogues of relations (|5.6[) and (|5.8[) . 

We next consider the general interacting case u ^ with finite n, addressing 
the issue of the series expansions of the Casimir amplitudes Ay m, n) and 

A^ BC (d, m, n) in e = d*(m) — d. We assert that these quantities have the Taylor 
and fractional power series expansions 

• ' ' ' ' ' C'(4) 5(n + 2) 



A 



1440 



and 



IE 



ln(2v^) - 



C(4) 4(n + ■ 



0(e 2 



(5.9) 



j-TT 2 


7r 2 e 


l \ 90 + 


180 


7r 2 (n- 




9-y/6 \n- 





IE 



ln7r 



2C ; (4) 5(n + 2) 
C(4) 2(n + 8) 



3/2 



,3/2 



+ o(e 3 / 2 ) 



(5.10) 



respectively. Here 7_e — ~r'(l) = 0.577215 ... is the Euler-Mascheroni constant, and 
C(s) means the Riemann zeta function. 

As is borne out by the second equation, the small-e expansion of A^ BC also 

involves half-integer powers of e. Moreover, beyond the given order e 3 / 2 , additional 
half-integer powers of e appear together with powers of lne. For example, there is a 
term proportional to e 5 / 2 lne. This is completely analogous to what was found in the 
CP case for A| BC [32l [33l [39] . The existence of such terms should become clear as we 
outline the derivation of the expansion (|5.10|) below. Before we turn to this matter, a 
few other remarks are appropriate. 

Relations (15.61) and (15.81) suggest to compare the above results with the analogous 



series expansions of A[,p' 0) and Ag BC in e = 4 - d. The comparison shows that all 
explicitly displayed orders in equations (|5.9p and (I5.10[) agree with their counterparts 
for A(^>'°^ and A™° obtained in [30j |3T] and [32j|33], respectively, up to the factor 
C m and the fact that e means 4 — d in the CP case rather than 4 + m/2 — d. 

An immediate important consequence of this finding and equation (|5.8j) should be 
mentioned. It was shown elsewhere [33] that the large-n limit of the fractional e = 4— d 
expansion of A^ c to 0(e 3 / 2 ) is in conformity with the e = 4 — d expansion of A^P^. 
Thus the large-n limit of the series expansion (|5.10l) in e = 4 + m/2 — d to 0(e 3 / 2 ) 
is also consistent with the small-e expansion of the large-n amplitude A B ^ (d, m) to 
this order in e. 



For example, the free energy of the Gaussian theory involves the function h(q 2 ) = lnq 2 
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Using relation (|5.8|) for PBC, we can determine A^(3, 1) in a straightforward 
manner by computing A^ c 00 (5/2). According to [10] and [33], ApP^ )0 (d) is given by 



a pbc 



(d) 



d Kd 



n,, 



Here Ad and the function Qd,2(y) are defined by 
A d = -(4n)- d ' 2 T(l-d/2), 

and 



(5.11) 



(5.12) 



E 



(d-1) 



(d) 
q=(p,k 



y 



d l 2 {2v\ 



d/2 



oo 

E 



-(d-2)/2 



K {d -2)/2{3Vy), (5-13) 



where K u (z) is the Macdonald function (modified Bcsscl function of the second kind) 
Further, Tvld is the solution to 

,(d-2)/2 



Adii^-*' 1 * = 21l- d l Q dt2 {n d ). 

Solving the latter equation at (i = 5/2 by numerical iteration yields 
TZ 5/2 = 0.5358 

from which the numerical values of the Casimir amplitudes 



A; 



PBC 



and 



cp,oo(5/2) = -0.2337. 
Afl^C 3 . 1 ) = -0.1269.. 



(5.14) 

(5.15) 

(5.16) 

(5.17) 
follow. 

It will suffice to give a brief outline of how the above small-e expansion results (|5.9[) 
and (|5.10j) were obtained. In order to determine Ay°'°' ) to 0(e), one can start from 
the two-loop expression for the total free energy F, set f = tlp and p = /3lp, subtract 
the bulk and surface contributions to F(ksTA)^ 1 to identify the residual free energy 
density f res (L; . . .), express the latter in terms of renormalized variables u and cr, 
and then evaluate f re s;R at the fixed-point value u* given in equation (|4.25j) and set 
L = a = /i = 1. If we proceeded in this manner in the case of PBC, we would 
encounter a zero-mode problem: The Fourier transform of the free propagator G£^ c 
at zero momentum (k,p) — becomes massless at the LP and hence IR singular. 
To see this, note that the spectral decompositions of the free propagators G^fjj at 
t = p = for the two BC in question read 



where 



(m) p(d—m—l) 



^ ^(z)[^(zT 



p 



2 + P 2 + ak 4 



exp 



d-l 

>E 

7=1 



(5.18) 



BC (z) 



(5.19) 



exp(iP r z), P r = 2nr/L, r <= Z, BC = PBC, 

y/2sm(P r z), P r = rir/L, r = l,...,oo, BC = (0,0), 

are the normalized eigenfunctions of the operator —d 2 on the interval [0, L], P 2 are 
its eigenvalues, and the sums ^ r run over the indicated BC-dependent sets. 
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For PBC, the lowest eigenvalue P% = Pq vanishes. Hence Landau theory predicts 
the film to have a LP at the same location (r, p) — (0,0) as the bulk. This is a 
deficiency of Landau theory. For (0,0) BC (corresponding to Dirichlet BC), this 
does not happen. Since the lowest eigenvalue Pp is strictly positive, none of the 
modes becomes critical at the bulk LP. 

To cope with this zero-mode problem, we proceed as in [32], [33] and [39]. We 
split off the (P r = 0) zero-mode contribution 

<p(y) =<5bc,pbc^ 1/2 / dz<fi(y,z) (5.20) 
Jo 

of the order parameter in the case of PBC, writing 

cf>(v,z) = L- 1 ' 2 V {y)+^{y 1 z). (5.21) 
Depending on the BC, one either has 

ip(x) = (j){x) and (p(y) = 0, BC = (O, O), (5.22) 



or 



/ dzip(y,z) = 0, BC = PBC. 
Jo 



(5.23) 



The field ip gives a contribution F^>, defined by 



-h^e-Wl, (5.24) 
k B T 

to the total free energy. We denote the associated reduced free-energy density as 

M L ^^Lkh- (5 ' 25) 

Its Feynman graph expansion to two-loop order is given by 

-U(L)A=( ) + W +0(3-loops), (5.26) 
* ' i i 

V / 

where the dashed blue line represents the free ^-propagator, i.e., the free propagator 
G^.y one obtains from the spectral decomposition (|5.18l) of Gf9, by omitting 
the (Pq = 0)-mode contribution in the case of PBC. Noting that the zero-mode 
contribution can be written in terms of free bulk propagator in d — 1 dimensions, 
one concludes that 

/•O) Ad-m-l) BCf \[ BC( ni» 

= Gl% (x; x') - fec.PBC L- 1 G { h d ' hm \y - y'). (5.27) 

The Feynman graphs displayed in equation (|5.26l) are computed in |Appcndix A| 

Since f-^(L) does not involve a zero mode, it is perturbatively well defined at the 
bulk LP. Its value at u = u* has an expansion in non-negative integer powers of e. 
For BC = (O, O), is identical to the total reduced free-energy area density f(L), 
whose definition should be clear by analogy with equation (|5.25[) . However, for PBC 
the contribution associated with ip must be added. We therefore write 

f(L) = U(L) + U(L) Wbc- (5.28) 



d-l 



i > g 7 (x 7 - x') 
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To determine f v (L), we substitute the decomposition (|5.21l) of <fi into the Hamiltonian 
and define an effective (d— l)-dimensional action H c s[(p] by tracing out ip. This yields 

c -u ctl[v ] = c F^/k B T Tr ^ e -«[*[* lV ]] (5.29) 

with 

y>]] = H[il>] +U[L- 1 ' 2 v\ + Hutlv, t/>], (5.30) 

where the contribution depending exclusively on tp and the interaction part are given 
by 

P P o m 2 1 ^ — 1 o 



E 



a=l 
2 "T 



/3=m+l 



a=l 



(5.31) 



and 



H in t [¥>, V] = ^ y [(^V 2 + 2( V • t/>) 2 + 2VT fo> • iA)^ 2 ] , (5.32) 

respectively. 

The effective action can be written in terms of an average (. . .)^ with the path 
probability density e F * _w ^2?[t/»] as 



In terms of the latter, the contribution f v {L) becomes 

U(L) = lim = - lim ilnTr^e-^^. 



(5.33) 



(5.34) 



Since we ultimately wish to use RG-improved perturbation theory, we need some 
results of the loop expansion of H e s Vp\ ■ 



n en [<p] = n[L-^ 2 (p] + H l £[<p] + 0(2-loops). 
The one-loop contribution is given by 

H [ £[<p] = I Tr In 



(5.35) 



W '1' Ifi 1 / \ / \ 

* W ^ U>' N - L 



<P. >- JP 

+ .... (5.36) 

Up to the order of our approximation (e 3 / 2 ) only the first graph explicitly shown in 
equation (|5.36[) contributes. It produces a shift f — > f + 5t^ c of the coefficient f in 

■H[L _1//2 y>] whose evaluation at r = p = and d — 4 + m/2 — e is straightforward, 
giving 



PBC 



n + 2 r(l- e/2)C(2- e) 



(5.37) 



6 27r 2 - £ / 2 L 2 - e ' 
The comparison with equation (14) of [32] shows that this shift is related to its CP 
analogue St^^, in the expected manner: 



Sr^ c (d, m) = a-™/4 Cm Jf PBC (d _ m/2) _ 



(5.38) 
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Both 0(u) shifts are proportional to £(2 — e), where £(s) is the Riemann zeta 
function. The latter is known to have a unique analytic continuation to the entire 
complex plane, excluding the point s = 1, at which it has a simple pole with residue 
1 [JTJ chapter 23]. This pole corresponds to a logarithmic IR singularity: When 
d — 3 + to/2 (e = 1), the free bulk propagator Gyf (x) decays at the LP for large \x\ 
so slowly that G^ B ^{x,x) diverges logarithmically. Thus it acquires an IR pole, just 
as the shifts (|5.37[) do. On the other hand, these shifts have an expansion about e = 
(i.e., about d*(m) = 4 + m/2 and eP(0) = 4, respectively). Upon expressing <Jt™ c in 
terms of renormalized variables and using £(2 — e) = 7r 2 /6 + 0(e), one arrives at the 
expansion 



r -PBC _ 
dT L,\\ - L 2 



7r 2 (n + 2) n . . 
— ^ L + 0(u,e) 



(5.39) 



Note that this quantity is positive when evaluated at the IR-stable root u* of 
P u . This is just as in the CP case: owing to this L-dependent shift provided by 
the coupling to the nonzero modes, the ip field does not become critical at the LP. 
Therefore, we can use the free propagator of the Hamiltonian (|5.31|) with p = and 
St^ c substituted for f in the Feynman graph expansion of f^ BC (L). This propagator 

is nothing else than the bulk propagator G^ 1 ' m \y\a; 5f, 0). 

Before doing this, a few remarks may be helpful to put things in perspective. 
The usual massless 4> A theory in 4 — e dimensions is well known to be IR singular in 
perturbation theory. IR poles occur at all rational e r = 2/r with r£N [44l l45l |46] , 
Nevertheless, the critical theory exists in the superrenormalizable case e > 0. The 
situation was clarified by Symanzik who showed that the bare critical mass parameter 
fcp (the CP analogue of tlp) is nonzero, nonanalytic in the coupling constant and 
not fully accessible by perturbation theory |441 145] . On dimensional grounds, it can 
be written as fcp = u 2 / £ T(e), where T(e) is a meromorphic function of e = 4 — d, 
with simple poles at e r . Likewise, one should have 

n jP =u 3 f e &- m ^T m (e), (5.40) 

where the functions T m {e) are expected to display analogous nonanalytic behavior in 
e = d*(m) — d at e r . 

In RG studies of CP properties based on the e expansion the problem of these 
IR poles is by-passed by choosing e smaller than the smallest value of the IR poles e r 
occurring at a given order of the calculation. Working at fixed dimensions d requires 
the use of massive RG scheme (see, e.g. 07], [IB], chapter 28.2] and the extension 
of this approach to systems with surfaces introduced by two of us [50] [51]). These 
schemes avoid the perturbative calculation of critical mass parameters such as top by 
expressing the quantities of interest in terms of the bulk correlation length. Crucial 
prerequisites for their quantitatively successful applications to the study of bulk critical 
behavior are the availability of perturbation-theory results to sufficiently high orders 
of the loop expansion [52[ [53] as well as detailed information about the large-order 
behavior of the required series (for details and references, see e.g. [49j and [54]). We 
are here in a much less favorable situation. First, we are not aware of any reliable 
quantitative investigations of bulk critical behavior at LP based on the massive field- 
theory approach, let alone of finite-size effects. (There is, however, some recent work 
using the so-called exact RG equations [SS].) Second, even the application of massive 
fixed-dimension RG schemes to the study of finite-size effects at CP has just begun 
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|56j . Third, the technical challenges one is faced with in such investigations of finite- 
size effects are enormous, even more so for systems with LP. 

Let us therefore take a pragmatic point of view and explore what the theory based 
on the smallness parameter e yields. To this end, we include the shift (I5.39[) in the 
free ^-propagator of the Hamiltonian (|5.31[) . The contributions to /„ we need to the 
order of our calculation then are simply given by the free-energy graphs to two-loop 
order of the (d — l)-dimensional action 

AU(L)= -Q - Q (5.41) 

where the red dot represents the coupling constant (—u/L). Upon exploiting again 
relation (15.51) , we can infer the resulting contribution to the residual free energy density 
fip,res(L) from equations (4.28) of [33]. We obtain 



Ad-m/2-l /^ PBC ^(2d-m-2)/4 



d-m/2- 1 



L. 



/j, e u n + 2 Ai-m/2-i /r.pBc^-Wa 

( dT L,\\ > 



(5.42) 



L 4! F ,, 

where Ad denotes the coefficient defined in equation (|5.12j) . and equation (|5.4p was 
used. 

The results given in equations (|5.39[) - (|5.42[) and (|A.7[) — (|A.8|) can now be combined 
in a straightforward fashion and substituted into equation (|5.28|) . Evaluating the 
resulting expression at the fixed-point value (|4.25[) of the coupling constant u then 
gives the small-e expansions (j5~9l) and (15T01) of A* ' ' and Ajp 30 , respectively. As 

can be seen from equation (|5.42p . beyond the order of e 3 ' 2 , higher half- integer powers 
appear in conjunction with powers of lne. 

In table [T] we present numerical values for the Casimir amplitudes Ay (d, m, n) 
of the Gaussian and n- vector models with d — 3 and m = 0,1,2. We do not consider 
(d = 3)-dimensional 4> A models with m = 3 because their lower critical dimension 
d*(m) = 2 + m/2 exceeds 3. However, the Gaussian critical case with m = d = 3, 
corresponding to a Hamiltonian with L\, = (A</>) 2 /2, is meaningful. Therefore, we 
have included it. 



5.2. Perpendicular slab orientation 

In our calculations of Casimir amplitudes for parallel slab orientations we could benefit 
from a simplifying feature: the required Feynman integrals could be related to their CP 
analogues. For perpendicular slab orientation, there is no such luxury. The necessary 
calculations become more involved, and especially in the case of (O, O) BC, require 
new techniques. 

We start by decomposing the total reduced free-energy area densities for PBC 
and (O, O) BC by analogy with equation (|5.28l) as 

fl C (L) = f™(L) + <5bc,pbc /*L(L) (5.43) 

into their nonzero-mode parts f^^j_(L) and a remaining zero-mode part f^ B ±{L) that 
is present only for BC = PBC. Using a subscript [I] to specify the contribution to 
f B( ~l(L) of Ith order in the loop expansion, we write 

f™(L) = f^ m {L) + /|£ ;[2] (£) + 0(3-loops), (5.44) 
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Table 1. Casimir amplitudes at d = 3 for parallel slab orientation. denote 
the Casimir amplitudes of the Gaussian models under boundary conditions BC. 
Results obtained for </> 4 models by extrapolation of the small-e expansions to d = 3 



are labelled by subscripts 'extr'. 


m 





1 


2 


A|^ c (3,m,n)/n 


-0.19131 


-0.15792 


-0.13090 


Af,^ tr (3,m,l) a 


-0.19670 


-0.14627 


-0.088776 


A||^; ^(3, m, n)/n 


-0.023914 


-0.027917 


-0.032725 


Ai° e 2(3,m,l)- 


-0.011659 


-0.0076387 


-0.0041162 




-0.15305 


-0.12698 




A ||,oo ( 3 ' TO ) 


-0.012(l) b 







a Values obtained by evaluating the small-e expansion at e — 1 + m/2 (d — 3). 
b Value taken from [38]. 



where the two terms on the right-hand side correspond to the analogues of the graphs 
displayed in equation (I5.26[) . 

We consider both BC, PBC and (O, O), separately, beginning with the technically 
simpler case of PBC. 

5.2.1. Periodic boundary conditions 

As before, we first consider the Gaussian case u = 0. Details of the calculation 
of the one-loop term m(L) are presented in |Appendix C| The result yields the 
residual free energy 

W - ° {d - m) ' 2 K B c7(d,m,n)L-^-™- D (5.45) 
with the Gaussian Casimir amplitude 

AP B C(ci,m,n) _ (0) 

Z _ Z7Tr rn,d+2 



(5.46) 



2 2d - 2m+1 T(d - m/2) C(2rf - m) 
~ 7r(d-i)/2r[(m-d+l)/2] ' 

where we introduced the quantity 

(0) = 2 2 ^ 2 "'- 4 r(d-2-m/2)C(2d-m-4) 

m,d ^-i)/2r[(m-d + 3)/2] 1 ' ' 

for subsequent use. In figure [2] the amplitudes (|5.46|) with m = 1,2,3 are plotted as 
functions of d. 

Turning to the small-e expansion, we note that the zero-mode contribution 
to /™^(L)|f. = ^ = o vanishes at zero- loop order. Hence the one- loop term 
f^res j_.fi] (^) j-7 L = y 5=o is also given by equation (|5.46[) . Both the two-loop contribution 

to an d the sn 'ft originating from the graph can be computed in 




a straightforward fashion along lines similar to those followed in |Appendix C| The 
results are 



T=p=0 



PBC 



n(n + 2) d _ 2 _ 



and 



Srf BC 



>L,± 



The quantity Pj^ d with m = 1 is plotted as a function of d in figure [3] Unlike the 
shift Srlf given in equation (15.39[). this quantity — and hence the associated St^ 3 ^ 
— is negative near the upper critical dimension d*(l) = 9/2. In the biaxial case m — 2, 

both Pjf?4 and the shift <5t^ j_ c vanish at the upper critical dimension d*(2) — 5. A 
nonpositive shift is a clear warning that a reorganization of RG-improved perturbation 
theory along the lines used in .32 , 33, 39] to deal with the zero-mode problem, and 
in section 15.11 in the case of parallel slab orientation, may not work in the present 
perpendicular case. Although one might hope that an appropriate extrapolation of 
the shift to d = 3 dimensions would yield a positive value, there is no guarantee that 
this is the case, nor that a meaningful small e-expansion will result. Unfortunately, 
we are not able to clarify these issues in a satisfactory fashion. We therefore take a 
pragmatic point of view, outline our strategy to derive a small-e expansion, and show 
that it fails when m > because the shift becomes negative for small e > 0. 



u a 



4! 



i /4 



3m/4 ^2m-4d+9 



P, 



(0) 



n(n 



m-2d+4 



P 



(0) 



(5.48) 



(5.49) 
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Figure 3. 



4.5 



5.5 



d 



3 (0) 



introduced in equation 115.471 1 and used throughout 



equations 115.481 1- 115.491 1. as a function of d for m = 1. 



The Hamiltonian - H[L~ 1/ ' 2 <£>] now differs from the expression given in (|5.31[) in 
that the sum over a is restricted to a = 1, . . . , m — 1 and the sum over /? runs over 
all d — m values /3 = m + 1, . . . , d. We make the replacement f — > f + 6t£^ in this 
Hamiltonian and set f = 0. Using again relation (|5 .5[) and remembering that the t/? 4 
coupling constant of this effective (d— l)-dimensional bulk action is given by u/L, we 
obtain 



PBC 
t£>,res;_l_ 



(L) 



, * -m/ 4+1/4 (~i 

1(7 O m _i 



Aj-ro/2-1/2 / A .pBC'\(2'i-m-l)/'l 

d- 7^/2-1/2^^ j 



(n+2) Ad-ro/2-1/2 /...pbcM 2 ^-™ 1 ^)/ 2 



4! 



■(ft 



L,_L 



(5.50) 



where C m and were defined in equations (|5.3|) and (I5.12[) , respectively. To compute 
the small-e expansion of the Casimir amplitude A^ BC , we must add f^^ ; ±{L) to 
the two-loop results for derived above, express the result in terms of 

renormalized quantities and evaluate it at u — u* , t = p = 0, and a = ji = 1. Since the 
shift c^t^ 1 ^ ~ e (in our attempted approximation), the first term in equation (15.501) 

varies as e 7 / 4 . Hence we expect a fractional power series expansion of the form 

A^ B (d,m,n) = + + , n y/i + f £ V*) (5.51) 

n 

whose higher-order terms should again be modified by powers of lne. 

The terms of zeroth and first order in e of this series originate entirely from 
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fK^.j_(L). The associated series coefficients a and a ± are well defined and given by 



29-m 7r (26-ro)/4 

an(m) ^- 1575r[( TO -6)/4] (5 - 52) 

and 

2 9-m 7r (26-m)/4 | n _ + In^^-S) + 3 ^,[( m _ 6 )/ 4 ] + 1 2 C'(8)/C(8)} 



ai(m, n) 



9450 r(m/4- 3/2) 

2 5-m 7] .(26-m)/4 ( m _ 2 ) 71 + 2 



(5.53) 



2025 r(m/4- 1/2) n + 8 
For the coefficient 07/4(771, n), one finds 

2- m + 7 7r 6 -"/ 4 r(-7/4) [(m-2)/4] 7 / 4 / / n + 2y /4 

a7 / 4(m ' n) = 3V 25 V4r[( m + i)/ 4 ] IrTTs J ■ (5 - 54) 

In the uniaxial case m = 1, it becomes complex, which renders the expansion useless. 
This is a direct consequence of the fact that St£± approaches a negative number as 
e — > 0+ when m = 1. 

5.2.2. Large-n results for perpendicular orientation and periodic boundary conditions 
As an alternative to the small-e expansion, we now wish to use the large-n limit to 
gain information about the Casimir amplitude A^ BC for the uniaxial case in d = 3 
dimensions. The result could be derived by solving a corresponding mean spherical 
model with short-range interactions along the lines followed in [40]. However, the 
large- n limit of the 4 model we worked with throughout this paper (cf. section [5]) can 
also be obtained directly. Since we are interested in the solution at the bulk LP, we 
set f and p to their bulk LP values tlp given in equation (|B.12[) and plp = and look 
for 0(n) symmetric solutions. Let us write the self-energy as — 1?^ ^ ■ In the large-n 
limit, one has i9 BB ^ = un(0 2 )/6 where the latter average must be computed with the 
self-consistent propagator involving the self-energy — i?^ B ^ (see e.g. [49]). To obtain a 
meaningful n — > 00 limit, we must scale ubyn, keeping 

g — nu (5.55) 

fixed. 

We can now subtract the self-consistent equation for and its analogue for 

the bulk quantity i?lp = ^lp° = li m L^oo at the LP to obtain a self-consistent 

equation for the difference 61}^ = i?^ B f - i?™ c . The bulk quantity i}^ c does 
not, of course, depend on whether the surface orientation is perpendicular or parallel. 
Furthermore, we can use the fact that the self-consistent bulk propagator is massless 
at the LP. We thus arrive at the equation 

r* p(m—l) p(d—m) 

9 ±,L 



1 ^ 1 



L p 2 + a [k 2 + (27rr/L) 2 ] 2 + ftjrao 
dK 1 



(5.56) 



, 2tt p 2 + a(k 2 + K 2 f 

whose analogue for parallel orientation, (|B.14I) , is used in |Appcndix B| 

The significance of the quantity <j$j^p(d, m; g) should be clear: it is the finite-size 
susceptibility, taken at the LP. We make it dimensionless, introducing 

H± = 6til B £L 4 /a, (5.57) 
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and then use Poisson's summation formula 

1 00 00 roc jtv- 

I £ /(2*r/£) = £ / (5.58) 

r— — 00 j ——oo 00 

to express the right-hand side of equation (|5.56D in terms of the free propagator G^ iW ^ . 
To this end, we define a function 

f (m— 1) p{d—m) 



Qi m (y) = § 



XSE27TZ " 2 " 



p 2 + (fc 2 + X 2 ) 2 + y 



4 d,m) 0e,|a = l; 2/^ = 0) (5.59) 



by analogy with equation (|5.13p . whose representation in the second line follows via 
the summation formula (|5.58[) . Evaluated at L = a = 1, the right-hand side of 
equation (|5.56|) then becomes 

2Qj; m (ftL) /- (m) /- (rf - m) gx , . 

A 7p (p 2 + fc 4 + 7^)(p 2 + fc 4 )- ^- bUj 

The dimensionless susceptibility 7£_l approaches a nontrivial number independent 
of L as L — » co. This number is the zero of the function (|5.60p . The term on the 
left-hand side of the self-consistency equation (15.56)) gives corrections to scaling, which 
we ignore. Upon evaluation of the bulk integral in expression (|5.60p we thus arrive at 
the self-consistency condition 

" Cm A d _ m/2 7^-W2-2)/2 + ^ Q ± m{n±) = o (5 61) 

for 1Z± . This applies to those values of d and m for which the required integrals and the 
function Qj[ rn are well-defined. This is the case if we require that 2 < d — m/2 < 4. 
For d — 3 dimensions, we are left with the interesting uniaxial case m = 1. The 
required bulk propagator 



(1) /■(<*-!) gjfcz 



p 2 + fc 4 + y 



' cos (zy 1/4 /^2) e-l z ^ 1/4 /^ (5.62) 



2tt|z 



can be calculated in a straightforward fashion. Inserting it into the last line of 
equation (|5.59j) and performing the summation gives 

QJ- ( y ) = _ JL i n h + e ~(4y) 1/4 _ 2 e -(«/ 4 ) 1/4 cos^) 1 / 4 ! . (5.63) 

47T L 

It is now an easy matter to determine the root lZ±(d,m) of equation (|5.6ip f° r 
(d,m) = (3, 1) by numerical means using Mathematica [57]. One obtains 

ft_L (3,1) = 0.998614.... (5.64) 

Thus the large-n limit yields a well-defined positive value for the scaled inverse finite- 
size susceptibility 1Z± at (d,m) = (3, 1), unlike the e expansion. 
To determine the Casimir amplitude A^ B ^ we must compute 

A F ± B £ = km £ 2d - m -7rcs,LP,ooU=i, (5.65) 

L—t-OQ 
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where the subscript oo again indicates that the limit / rC s,LP,oo = ri m n->-oo /res lp/^ 
has been taken. This limit is given by the value D(JZ±) of a function D(y) that 
satisfies D(0) = A± B £{d, m, 1) and can be determined by integrating the equation of 
state (|5.6ip with respect to 1Z± . Hence we have 



A 



PBC 



(d, m) 



A 



d-m/2 „ „(2d-m)/4 



aPBC 



(d,m, 1) 



^Qi m {y). (5.66) 



u-m/2 Jo 

We now set (d,m) = (3, 1), insert Q^iiu) from equation (I5.63[) . and use Mathema- 
tica [57] to perform the integral. This gives, after some rearrangements, 

5 / 4 V2 



ATO3,1) 



_ — Re^ (1 



Siv^y 1 /2 Li 3 



*)y 3/4 Li 2 



(6 + 6i)y 1/4 Li 4 



6V2L 



15 



(l+ l )(y/4) 1 / 4 



!/=Ki(3,l) 



~ 4.00053, 

where hi s (z) is the polylogarithm, defined by [5 8) 



Li s (z) 



OO 

E 



(5.67) 



(5.68) 



Note that the sign of this Casimir amplitude is positive. That is, the Casimir force 
is repulsive in this case. This is remarkable because, under PBC, the corresponding 
Casimir forces at uniaxial LPs and parallel orientation and at CP are both attractive 
[cf. equation (|5.17j) and table [IJ on the one hand, and equation (|5 . 16(1 . on the other 
hand] . More generally, there exist theorems [42l [43] which state that Casimir forces 
are guaranteed to be negative (attractive) when certain symmetry properties (such 
as the same BC on both planes and reflection-positive interactions) are satisfied. In 
the LP case with perpendicular orientation, reflection positivity is violated along the 
z-axis due to the presence of competing nearest and axial next-nearest neighbour 
interactions along this direction. Thus the theorem of [43] does not apply. The one of 
[42] is restricted to Gaussian models. It also does not apply to our LP case, neither 
on the level of the corresponding Gaussian model nor in the large-n limit, because it 
presumes a Gaussian action involving square-gradient terms of the order parameter 
but no higher-order derivatives. 



5.2.3. Ordinary- ordinary boundary conditions 

The case of perpendicular orientation with (O, O) BC turns out to be technically 
very demanding. For the sake of simplicity, we will therefore content ourselves here 
with an analysis to one-loop order. 

To compute /„^'°?m (L), we proceed as follows. We consider the Gaussian infinite- 
space model with Hamiltonian 

H G = [ d d z£ b (aO|. =0 .. . . (5.69) 

We now insert two defect planes into the system, one at z — and a second at z = L, 
where z is taken to be an a-direction. At either one of them the BC (|3.5[) are required 
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to hold. We then compute the change in reduced free energy AFq^bT) 1 caused by 
the defect planes. The L-dependent part of the area density lim^oo Fq^bTA)^ 1 
gives us the desired residual reduced free-energy density /^'^j (£). Note that this 
conclusion hinges on the BC p.5[) , which ensures that the region between the defect 
planes is decoupled from the outside regions z < and z > i|j 

We implement the BC (13.51) at the planes *8i (z = 0) and < B 2 (z = L) in the 
required functional integral by delta functions. Using a standard trick |61[ 162 ) . we 
represent these constraints through functional integrals over four n-component fields 
•dj whose supports are restricted to *Bi and < B 2 for j = 1,2 and j = 3, 4, respectively. 
Thus 

^(xiM^M J] *[0(*2)]«[0„0(* 2 )]= [v[#]c- [C ^ (5.70) 
xis»i x 2 e < s 2 

with 

C[0,tf]= / [-& 1 -<t> + -d 2 -d z <j>]+ ( [d 3 -<f> + 4 -d z <p], (5.71) 

where the functional integration J 2?[i9] is over all four fields 
As a straightforward consequence we obtain 

p -A W= Z't W 2 ^ exp(-ft[<fl-iC[«M]) 
7 U jVWexpi-Hfr]) 



T>[&] exp 



^<(C[0^]) 2 )« C 



(5.72) 



where (. . -)h g denotes an average over cf> with the path probabililty density cx 
e~^ G T> [<fi]. The remaining Gaussian functional integration over •& is straightforward, 
giving 

/ G (L) = lim 4i% = S (const + lim A' 1 In det Ml) , (5.73) 

A— >oo AKb± i \ A— >ca / 

where "const" is a constant which will drop out in the residual free energy 

f G ,rUL) =/ g (£)-/g(oo) (5.74) 
we are interested in. The matrix Ad l can be written in a suggestive manner as 
( i\G%' m \ i|G^ m) ^li ilG^b i\G { h d ' m) t z \ 2 \ 



Ml = 



alG^li 2|G^ m) K|i 2 |G^ m) | 2 2 \G^ m) t z \ 2 
V 2|^,G^ Tn) | 1 2 |^G^ m) K|i 2\t z G^ m) \ 2 2 z G { h d ' m)< d z \ 2 J 



(5.75) 

It consists of blocks involving the free bulk propagator G^'"^ and its derivatives 
at points on the defect planes Bj. The derivatives d z and a z act to the left and 

* In general, systems with two defect planes must be carefully distinguished from slabs confined by 
two parallel plates at a distance L. A decoupling analogous to the one just mentioned occurs in the 
case of the standard infinite-space </> 4 model in the presence of two defect planes at which Dirichlet 
BC hold 1351 . Dirichlet BC would also suffice to achieve this decoupling in the LP case if the defect 
planes were oriented perpendicular to a /3-direction. For a detailed discussion of the similarities and 
differences of boundary critical behavior at free surfaces planes and defect planes, the reader might 
want to consult [59], [BO] and [28] chapter IV.A.4]. 
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right, respectively. Further, the notations j\ and \ji indicate that the left point 
x i = (U2, z i) lies on *Bj and the right point x 2 — (2/2^2) is located on 23j/. 
For example, the elements of the block i\o z Glf' m ^ d\ 2 (labelled by y\ and y 2 ) are 
9zi d Z2 G^" m) [(yi ,zi=0),(y 2 ,z 2 =L)}. 

To exploit translation invariance along the y-direction, we Fourier transform with 
respect to y. Let (k,p) with k £ K m_1 and p € R d_m be the wave-vector conjugate 
to y. In this momentum representation, the matrix M-l is block diagonal. The 
logarithm of its determinant becomes 

p[m— 1) p(d—rn) 

lim A- 1 hxdet Ml = / / lndet A4 L (fe,p), (5.76) 

where .Ml(A;,p) is a 4 x 4 matrix. To compute its elements, it is helpful to know 

Md,m) n s f 00 di^ exp (ii4T(zi - z 2 )) 

k (fc ' P;Zl - Z2)= i_ 00 27 ^ + a(P + ^ ' (5J7) 

the free bulk propagator in the fcpz-representation. It may be gleaned from equations 
(20)-(23) of [26]; we have 

,(i,m) a m ^ _ k- cos(K_g~ 1 / 4 |z|) + k + sin(K_o-~ 1 / 4 |z 

4a 1 /4 K _ K+ ( K 2 +K 2_) 

with 

1 

V2 

A straightforward calculation then yields 

~ k 2 _ cosh(2£a- 1 / 4 K+ ) + 4cos(2L ( 7- 1 / 4 K _) - n 2 _ - K \ 

aetML\k,p) = — „ — j-r^ tt^ ttttz — rn — ; ■ (5.80) 

y ' 128<t 2 ^k 4 (k 2 _ + n 2 + ) 2 exp(2La- 1 / i K+) v ; 

Noting that 

lim detM L (k,p) = [256<7 2 k 4 (<L + /si ) 2 ] _1 , (5.81) 
we can make the subtraction in equation (|5.74j) to obtain 

/G,res(-^) 



mJ (fc,p;z) = K A «_ UA r U J, _2 .2 ^e~^- L '^ (5.78) 



= - v vV + & k 4 T o" 1/2 fc 2 - (5.79) 



(d— m) p(m— 1) 

In 



«£. cosh(2ff- 1 / 4 LK + ) + k 2 cos(2(7- 1 / 4 La 



-2 1 



(rei/2) exp(2<7- 1 /4i K+ ) 



(5.82) 

Assuming that m > 1 and d — m > 0, wc perform the angular integrations. This leaves 
us with the two radial integrals over p and k. To simplify them, we make the change 
of variables 

fc = 7&' p = ^' dfcdp = ^pr d ^ ( 5 - 83 ) 

where 



We obtain 



w, : \l x 2 + \Jx± + 1 . (5.84) 
/ Gj r es (i) = a^-^' 2 Agi 0) (d, m, n) L 1+m - 2d (5.85) 
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r ,2d-m-1 



with the Casimir amplitudes 

TS ts />oo „jnn— 2 roc 

A O.O), , W K d _ m K m _ 1 A X A 2 

A^'(d,m,n)/n= / d^ _ x / dgg 

z ^ " Jo U>x Jo 

x In (2e- 9 {cosh(g) -l+wi[ cos(g/u£) - l] }) . (5.86) 

Since <r = tr when u = 0, we have replaced <r by a in equation (|5 .85|) . 

This is our final result for A^P^ '' when m > I and d > m. In the special cases 
m = 1 (uniaxial LP) and d = m (isotropic LP), the result simplifies to single integrals; 
one finds 



A { °f\d,l,n)/n = Ka-^ 1 -* / dqq 2d ~ 3 In {2e~ 9 [cosh(q) + cos(q) - 2]} (5.87) 

Jo 

and 

poo 

A i °f\d,d,n)/n = 2- d K d - 1 dq q^ 2 In {e~ 9 [2 cosh(g) - 2 - q 2 ] } . (5.88) 

Jo 

The simplest way to get these results is by going back to the appropriate analogues 
of equation (|5.82l) . To obtain the analogue for the uniaxial case from this equation, 



- (m— 1) 



one must drop the integral J^ m and replace the integrand by its k — > limit; 



likewise, one must omit and take the limit p^Oof the integrand in the case 

of the isotropic LP. The result (|5.87[) for the uniaxial case can also be inferred directly 
from the general one, equation (|5.86p . with the aid of the formula (see e.g. [2"T] ) 

f°o r( D ) 
MmK D &KX D - 1 f D {x)=\mx I /d(x) = / (0), (5.89) 

where fjj is a rotationally invariant function of x £ H D that depends parametrically 
on D. This works because the original limit k — > at fixed p translates into the limit 
h 0. Deriving the result (|5.88[) for the isotropic LP from equation f|5 . 86[) is possible, 
but more subtle. (One would have to take the appropriate scaling limit q ~ p and 
h ~ p~ x / 2 corresponding to the original lim p „>.o of the integrand.) 

In table [2] we assembled various values of A.\!?±\d, m,n)/n for different m at the 
corresponding upper critical dimension d* = 4 + m/2 and for general m < d < 8. For 
the case of d = 3, we included all values for m < 3, in particular the special ones for 
the uniaxial (m = I) and the isotropic (m = 3) LP. 



Table 2. Casimir amplitudes Aq ' ± ' for various values of d and m. 
m 12 3 4 

A[?; 0) (d*,m,n)/n -6.596 -2.841 -1.271 -0.5635 

A { °f\d*,m + 4,n)/n -0.2411 -0.09856 -0.03848 -0.01438 
A^°i 0) (3,m,n)/n -1.3240 -0.5900 -0.2732 



6. Summary and conclusions 

In this paper we studied the effects of confining anisotropic scale-invariant systems 
between two parallel planes at distance L. Just as in the case of confined isotropic 
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scale- invariant systems, the confinement of long wave-length fluctuations induces long- 
ranged effective forces. However, in the anisotropic case, the problem is much richer. 
Important qualitative and quantitative differences arise. 

A first important qualitative difference implied by the distinct scaling behavior 
of coordinate separations along a- and /3-directions is that the orientation of the 
boundary planes matters. The exponents £j_ and £y characterizing the algebraic 
decay of the fluctuation-induced interactions between the boundary planes as a 
function of L differ, depending on whether the surface planes are aligned perpendicular 
to an a-direction or parallel to all a-directions, i.e. perpendicular to a /3-axis [cf. 
equation (jl.4|) ]. 

A second difference is that the proportionality constants appearing in the 
analogues of the decay law (|1.3|) involve nonuniversal amplitudes. In order to obtain 
universal Casimir amplitudes, one must split off such nonuniversal factors. As 
expounded in section l4~5I this can be achieved in a natural fashion by defining universal 
Casimir amplitudes as universal amplitude ratiosQ For given perpendicular 

or parallel surface plane orientation, these Casimir amplitudes depend on gross 
surface properties such as boundary conditions (BC). In cases of free BC, the BC 
that hold asymptotically in the large-length-scale limit are associated with the RG 
fixed points governing the surface universality classes of the respective semi-infinite 
systems. Owing to the different scaling behaviors of distances along the a- and /in- 
directions, the Hamiltonians of appropriate minimal continuum models involve distinct 
boundary contributions £j~3y 2 m the cases of perpendicular and parallel surface plane 
orientations. We introduced these models and discussed the associated mesoscopic 
BC in section [21 

Unfortunately, full analyses of these models for general values of the surface 
interaction constants are not even available for semi- infinite systems f22j [23l [24J [25] . 
For the sake of keeping the technical difficulties manageable, we restricted ourselves 
to the study of asymptotic large-length-scale BC of the (O, O) type when considering 
free BC, namely the Dirichlet BC Q3.2p and the BC <j> = d n <p = [equation p.5[l ] 
for parallel and perpendicular surface-plane orientations, respectively. In addition, 
we investigated the cases of periodic BC (PBC), both for parallel and perpendicular 
surface-plane orientations. The technical difficulties we encountered turned out to be 
somewhat easier to handle for parallel surface-plane orientations. For this orientation, 
we were able to compute the Casimir amplitudes Ajj^°^ (m, d, n) and An > Q C (m,d, n) 
of the corresponding Gaussian models. The results are given in equation (|5.2p . 
According to equation (|5.6j) . these quantities are related in a simple manner to their 
CP analogues A||(A re) and A^§ G G (D,n) at D = d - to/2 dimensions. As we 
showed in |Appendix B[ the analogous relationship, equation (|5.8[> . holds between the 
large-n Casimir amplitudes A^^ o (d, to) and Aj?p ^(d — m/2) for (O, O) and periodic 
BC. Our numerically determined exact value for AP^ (3, 1) is given in equation (|5.17j) . 

We also investigated the Casimir amplitudes Aff > '°\d,m,n) and Ay >BC (cZ, m, n) 
by means of RG- improved perturbation theory in d = 4 + to/2 — e dimensions. As 

(I The issue of universality was recently considered in some detail also for the Casimir amplitudes 
of weakly anisotropic scale-invariant systems [561 119) . As can be seen from these references (see, in 
particular, 1191 p. 15-17]), the weak anisotropy can be absorbed by a proper choice of variables. Such 
a transformation to a system with a standard square-gradient term in the action's bulk density will 
in general change the boundary interaction constants of the corresponding continuum model and, in 
cases of fully finite systems, also their shapes. 
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expected from the analogy with the CP case, , m, n) turned out to have 

a Taylor series expansion in e. To order e, it is given in equation (|5.9p . Owing 
to the presence of a zero mode at zero-loop order, the conventional RG-improved 
perturbation theory for A FBC (d, m, n) becomes ill-defined at the bulk LP because of 
IR singularities. To cope with this problem, we proceeded along the lines followed in 
the CP case [32l |33l [39] and in analogous problems of finite-temperature crossovers 
near quantum critical points [63] , using RG-improved perturbation theory to construct 
an effective (d— l)-dimensional action for the zero-mode component ip(y) of the order- 
parameter field (p(y,z). The resulting modified RG-improved perturbation theory 
showed that the small-e expansion of A PBC (d, m, n) becomes a fractional one of a form 
similar to that of its CP analogue A^p (d,n). It involves, besides integer powers of 
e, also half-integers powers t k / 2 with k = 3, 5 . . . that are modulated by powers of In e 
when k > 3. 

Performing analogous calculations for perpendicular orientation turned out to be 
technically considerably more demanding. For the sake of simplicity, we therefore 
restricted ourselves to a one-loop approximation when considering the case of the 
(0,0) BC (|3.5j) . This yielded the exact analytical expressions (|5.86[) . ()5.87[) . and 
(|5 . 88[) for the Gaussian Casimir amplitudes Aq 1 ^ (d, m, n) in terms of double or 
single integrals. In the case of PBC, we were able to perform a two-loop calculation, 
but were faced again with the problem that a zero mode is present in Landau theory 
for finite L at the bulk LP. To deal with it, we pursued the above-mentioned strategy 
of constructing an effective (d— l)-dimensional action for the zero-mode component of 
the order parameter by integrating out the remaining nonzero- mode degrees of freedom 
by means of RG-improved perturbation theory. Unfortunately, the mass term of the 
free propagator of the resulting effective action turned out to become negative near the 
upper critical dimension in the physically interesting uniaxial case. Upon inspection 
of the e dependence of the corresponding formal perturbation series of A PBC (d, m, n) 
we could identify terms involving fractional powers ~ e 7 / 4 and analogous fractional 
powers of higher orders times powers of In e [see equations (I5.51[) - (|5.54[) ] . Owing to the 
mentioned behavior of the free propagator's mass term, we refrained from attempts 
of extrapolating this series to d = 3. (Such extrapolations would require reasonable 
assumptions about the resummation of the series.) On the other hand, we found 
that a well defined large-n limit of A™ c (3, l,n)/n exists; the exact limiting value 
Ajj^j(3, 1), which we determined by numerical means, is given in equation (I5.67[) . 
Interestingly, this quantity is positive, so that the corresponding Casimir force is 
repulsive. 

To put our results and the problems encountered in the applied RG-improved 
perturbation theory in perspective, it is appropriate to recall a general difficulty one 
is faced with in studies of films near bulk criticality in the limit of large film thickness L: 
ultimately, one would like to have a theory that can handle the involved crossovers from 
(i-dimensional (multi) critical behavior to (d— l)-dimensional (multi)critical or pseudo- 
(multi) critical behavior. However, perturbative RG approaches that are capable of 
dealing with critical behavior in bulk and semi-infinite systems (such as those based 
on the e expansion) fall short of reaching this ambitious goal. Even the more modest 
aim of obtaining well-defined asymptotic expansions in e at the bulk (multi) critical 
point may prove elusive for certain large-distance BC involving zero modes at the level 
of Landau theory. It comes as no surprise that these problems turned out to be more 
pronounced and severe than in previous analogous studies of Casimir interactions at 
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bulk critical points [301 1311 1321 1331 13S] • The reason is that the present study of finite- 
size properties involves a bulk multi critical point, the LP. If the values of d and m are 
such that a LP continues to exist for finite L, the theory must be able to account in a 
reliable fashion for the shifts of the LP and the phase boundaries that meet at it. For 
d = 3 and finite L, the situation is even more complicated. One expects the thermal 
fluctuations to destabilize the phases with homogeneous and modulated order so that 
the critical lines get replaced by pseudo-critical ones and no real LP may be expected 
to exist for finite L (see [12j[l3j|64] and their references). 

An appealing feature of the large-n limit is that it is capable of dealing with 
thermal fluctuations and dimensional crossover in a mathematically controlled fashion. 
On the other hand, one must keep in mind that in the large-n analyses presented 
here we fixed the thermodynamic fields to their values at the bulk LP point and 
hence did not use them to investigate the phase diagram for finite L. The results 
are meaningful for the particular (large-n or spherical) models to which they apply. 
However, caution is necessary in applying them to the case of finite n. To draw 
meaningful conclusions, knowledge about the phase diagram from other sources such 
as Monte Carlo simulations or alternative approaches is necessary. 

Finally, let us briefly comment on the question of verification of the Casimir 
amplitudes investigated here. Previous experimental investigations of thermodynamic 
Casimir forces have focused for very good reasons on fluid systems [5J [5J ESI EMI ES] ■ 
The advantage of using fluid systems is that the number of degrees of freedom 
of the medium (fluid) between macroscopic bodies (e.g., walls) can vary, a fact 
which enables direct and indirect confirmations of such fluctuation-induced forces. 
By contrast, the continuum models exhibiting anisotropic scale invariance on large 
distance scales considered here have natural realizations as lattice models (such as 
the ANNNI model in the uniaxial case [12j [T3J [64]) and rely on the presence of 
corresponding lattice anisotropies. Unlike fluid systems, such systems do not lend 
themselves to direct measurements of thermodynamic Casimir forces. On the other 
hand, the Casimir amplitudes we considered are finite-size quantities. Just as other 
finite-size quantities (which are not necessarily measurable via induced forces) , they 
are well-defined observables that at least in principle should be measurable. At the 
moment, the most promising means of checking our predictions appears to be Monte 
Carlo simulations. This technique was successfully used already some time ago to 
study the surface critical behavior of bounded three-dimensional ANNNI models 
[68] . Appropriate extensions of this work along the lines of recent Monte Carlo 
investigations of thermodynamic Casimir forces near critical points [5H1 [701 [H] should 
be possible and enable reliable checks of our predictions. We hope that the present 
work will stimulate such investigations of fluctuation-induced interactions in strongly 
anisotropic critical systems. 
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Appendix A. Feynman graphs for the case of parallel slab orientation 



In this appendix we compute the one-loop and two-loop free-energy Feynman graphs 
that are needed to obtain our results for parallel slab orientation presented in 
section 15.11 

To determine the one-loop graph in equation (|5.26l) . we set f = p = 0, evaluate 
the required trace in the eigenfunction representation, and then exploit relation (|5.5I) 
to get rid of the quartic fc-dependence. This gives 

5 2 H[ip] 



-1 



J /A 2A Tl ' ln [6ip a (x)6Tp a/ (x>) 



il>=0. 



—n 
~2~ 



(m) p(d—m—l) 
k J p 



pi 



-n 



-m/4 s~i 



ak 4 ) 

Hi?]: 



(A.1) 



£ln( P 2 

P r ^0 

ro/2-1) oo 

where P r are the BC-dependent momenta specified in equation (|5.19l) . The prime on 
Tr reminds us that the trace in the case of PBC is restricted to the subspace orthogonal 
to the P r = mode. 

The integrals over the series in the last line are known from [31 j . |33j . and |35j . 
Using these results one arrives at 



AgS G (d 



nLX h (d - m/2) + nX® c (d - m/2) 



j/2) L-id-l-m/2) 



(A.2) 



where A^p G (D) is the CP Casimir amplitude of equation (|5.6p at dimension D = 

d— m/2. Further, X\>{D) {— J^ D ln/c 2 ) and X® C (D) are L-independent functions 
of D that vanish in dimensional regularization but would diverge as A D and A 13-1 , 
respectively, if a large-momentum cutoff A were used to regularize the UV singularities. 
Both terms do not contribute to the residual free energy density. The L-dependent 
last term yields the contribution to / r ^g (X; 0, 0, 0, a, fj,) given in equations (I5.ip and 

The two-loop term in equation (|5.26[) reads 



,5« j A 



n(n + 2) u 



(A.3) 



Upon substituting the eigenfunction representation of the propagator G^^.y, shown 
in the first line of equation (15 -27[) . we can once more harness relation (|5.5[) to conclude 
that this graph is given by its CP analogue in D = d — m/2 dimension times the 
factor a~ m / 2 C^ n . Insertion of the results for the corresponding CP graphs given in 
the literature [32 133] then yields 



i i 

where lf c (D) stands for 



n(n + 2) u C 2 



-h-2d+m 7-BC 



8 a m / 2 ' 



I^{d-m/2), 



li BC {D) = \2- x -k- d I 2 C,{D - 2)r(£>/2 - 1) 



(A.4) 



(A.5) 
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or 

AO,Q) (m g^V^[2 C 2 (3 -D) + C(6 - D)\ 
2 { >- T*[{D-1)/2]cos>{Dk/2) ' [ b) 

depending on the BC. 

Insertion of these results into equation (|5.26p gives 

= £ji ^lloid - m/2, n) L"*Hn/W 

+ n(n + 2) I ^_ 7 BC (d _ m/2) L -2 d+ro+5j BC = pBCi (Qj Q) ; 

(A.7) 



where A^p G (D,n) are the Gaussian CP Casimir amplitudes 



Appendix B. Large-n limit 



The purpose of this appendix is to prove that the relation (|5 . 6[) between the Gaussian 
LP Casimir amplitudes A^q and their CP analogues A^p G carries over to the ti^O 
theory in the limit n — > oo, as stated in equation (15.81) . 

To study the models introduced in section [2] in the large- n limit, we rescale the 
coupling constant as in (|5.55l) and consider the partition function 

ii;] S Jv[4>]e-«M. (B.l) 

Using the Einstein-like conventions that pairs of indices a, a and ft are to be summed 
over the values a = 1, . . . , n, a = 1, . . . , m, and ft = m + 1, . . . , d, respectively, we 
write the Hamiltonian as 



Z n = exp 



where K, means the differential operator 



^ 4 
4!n 



/c = *(di) 2 



(B.2) 



(B.3) 



subject to the boundary conditions discussed in section [3] We do not have to specify 
these at this stage but will do so later. 

Following a standard technique |72j . we introduce the energy density operator 
£ (x) = 4> 2 {x) and rewrite the interacting part at a given position x as 

9 „ 



exp 



/oo 
d£5(n£ - (j) a (j) a /2) exp 
-oo 



gn £2 



6 



n 

2xi 



di) 



d£ exp 



E 2 



(B.4) 

We insert this expression into the partition function (|B.1|) and then do the functional 
integral over <f> first. This is Gaussian, giving 



J 2?[0]exp j-i^[0 a (/C 



exp 



lndet (K + f + 4 - pd 2 a ) + - N x \n(2w) 



(B.5) 
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Here N x — v Q 1 J m dV is the number of lattice points x, where vo denotes an 
elementary unit cell (discretization volume) . The partition function becomes 



C n / X>[0] / V[£] exp 



1 



6 f 2 J--lndct (£- 



■ (B.< 



with C n = {n/\) Na > (27r)( n ~ 2 > N <°/ 2 . To obtain the large-n limit, we can evaluate the 
functional integrals by the method of steepest descent. The saddle-point equations 
become 



§{x) 



fax), 
1 



E{x) = -(x\{K + T + d-pdl)- l \x). 



(B.7) 
(B.8) 



Let £ (x) and J?(x) be solutions to these equations that maximize the integrand in 
equation (|B.6[) . Then the large- n limit of the reduced free energy per number n of 
components is given by 



lim 



F„ 



1 



IndetlX + i + tf-jod 2 ) 



£- 



ln(27r) 
2«o 



.(B.9) 



For bulk systems and finite-size systems with PBC, the saddle-point 
equations (|B.7J) have spatially homogeneous solutions (%9(x),£(x)) — (i?b>£t>) and 
($P BC , ££ BC ), respectively. Furthermore, r^ — i?t + f and = +r are the inverse 



bulk and finite-size susceptibilities. At the LP, both rb and the k = q a term of the 
Fourier transform of the bulk two-point vertex function K, + rb — f> d\ must vanish. 
Hence the critical values are given by 

(B.10) 
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If we continue to use dimensional regularization (assuming that d>2 + m/2to avoid 
IR singularities), we can use again equation (|5.5[) to relate tlp to its CP analogue in 
D = d — to/ 2 dimensions: 

f LP (d, m) = a"™/ 4 C m fcp(d - m/2), 
q 1 

tcp(D) =- 9 - - (B.12) 

6 J q r 

We now set (t,/3) = (tlp,Plp)j introduce the deviations of d and £ from their 
bulk LP values, 

8&{x) = $(x)-<d hV , 8£(x)=£(x)-£-l F , (B.13) 
and consider the case of parallel orientation for PBC. Looking for spatially 
homogeneous solutions S l 3^^(d, m; g) and 5£ jf B (d, m; g), we find the self-consistent 
equations 
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We can now apply relation (|5.5I) to the momentum integrals J^ d of the 



function h(p 2 + a k ) inside the square brackets. The resulting transformed equation 
agrees with its counterpart for the CP quantity 8 l d^ c L {D;g') of a (D = d — Tri- 



dimensional system with the coupling constant g 
have 



S^ B L c (d,m;g) 



J CP,L 



ga- m ^ C„ 
(d-m/2;ga- m / 4 C m ) 



Consequently, we 
(B.15) 



and 



0t \\,L 



: (d,m;g) = a^Cn 6£™ c L (d - to/2; g&- m ^C m ). (B.16) 
In view of equations (jB.111) and (IB.12|) , it is clear that these relations carry over 



to ^|"^ (d,m;g) and £ F B (d, m; g). Hence the term oc £ 2 in equation (|B.9[) can be 
expressed in terms of its CP analogue as 



£jT{d,m-g) 



° Lm ~Q V t -CP,L\ D -i9 )\ 
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That the Gaussian free energy term in equation (|B.9|) is related to its CP analogue 
in the same manner follows from our analysis of the Gaussian theory in section 15.11 in 
conjunction with the results of | Appendix A| and relation (|B.15[) . The upshot is that 
the n = oo Casimir amplitudes Ay'^ are indeed related to their CP analogues A™ 1 ^ 
as stated in equation (|5.8[) . This conclusion exploits the fact that the dependence on 
the coupling constants g and g' drops out of these Casimir amplitudes because of their 
universality, albeit giving rise to usual corrections to scaling. 

Next, we turn to the case of parallel orientation and (O, O) BC. Owing to the 
breakdown of translation invariance along the z-direction, the corresponding solutions 
$11^°^ and to the saddle-point equations (|B.7j) - (|B.8|) become z-dependent. The 

equations for the shifted quantities 6&\?j°^ and S£ n^' ' now read 



8^°l°\z\d 7 m-g) = | <5£|j°' 0) (z\d, to; g) 



... L 

(m) 



k 



(d-m-l) 



E 



\fr,L(z; g)f 



£r,L(9) 



k 4 



dP 



1 



2tt p 2 + P 2 + a k 4 



(B.18) 



Here {tp r ,L,} is an orthonormal complete set of eigenfunctions on [0,L] (which can be 
chosen real) and e r ,L are the corresponding eigenvalues. These functions satisfy 



-d 2 z +o\° L °\z;g) 



£r,L 



¥>r,L(z;g) = o 



(B.19) 



and Dirichlet BC at z 
read 



and z = L. Their completeness and orthonormality relations 
^ < L (z; g) cp r , L (z'; g) = S(z - z') (B.20) 



and 



respectively. 
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Obtaining analytical solutions to the above self-consistent equations is quite a 
challenge. Whether this goal can be achieved is unclear to us, both in the LP 
case we are concerned with here and the CP case with Dirichlet (or free) boundary 
conditions. Bray and Moore [37] managed to find analytical solutions for the CP 
quantity (5i9cp.l=oo(z) in semi-infinite geometry. However, for finite L, one generally 
will have to resort to numerical techniques, as was done in a recent large-n study of 
the CP Casimir effect [38] . 

Numerical calculations of the eigensystem {e r ,L, fr,L} require in one way or 
another a cutoff A z = 7r/a z , where a z is a discretization length (lattice constant) 
for z. On the other hand, it is preferable to do the required momentum integrations 
j(d m i) j{m) ana jy^j ca Uy as before, and not restrict them by a cutoff. The price is 
tnat UV singularities are encountered in the free energy, even for d strictly below the 
upper critical dimension d*(m) = 4 + m/2 (where the theory is superrenormalizable) . 
To obtain convergent integrals, one must (and can) make appropriate subtractions. 
(Subtracting the bulk and surface free energy contributions to obtain the excess free 
energy density will remove the bulk and surface UV singularities associated with 
additive counterterms. For d < d*(m), the resulting differences must be UV finite.) 
In the equations for the shifted variables, UV fmiteness is ensured when d < d*(m) 
because of the subtracted bulk terms. Reasoning as in the case of PBC, we can 
therefore conclude that the following analogues of relations (|B.15|) and (|B.I6|) hold: 

5^ L °\z\d,m-g) = 6^°^(z\d-m/2;ga- m / 4 C m ) (B.22) 

and 

8£l°<°\z\d,m;g) = ^6e&>°\z\d-m/2-,g&- m ' i C m ). (B.23) 

The residual free-energy densities fft®p (L\d,m; g) and f^p'°\L\D- g') must 
therefore be related in the same manner, which in turn implies relation (|5-8[) for 
the Casimir amplitudes. 

Appendix C. Feynman graphs for the case of perpendicular slab 
orientation and periodic boundary conditions 

In this appendix we compute the Feynman graphs required for our calculation of the 
Casimir amplitudes for perpendicular orientation and PBC reported in section f5. 2 .11 
We begin with the calculation of the one-loop contribution to the free-energy area 
density 

f(d-m) Am-l) 

flf ] (L\a;f,p) = - / $>fi,(g P |f,p) 
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for f = p = 0. Here 

to*{qr\T,f>) =T+p 2 +pk 2 + cr[k 2 + (2irr/L) 2 ] 2 (C.2) 

and q r denotes the d-dimensional momentum vector (fci, . . . , fc m _i, • ■ • ,Pd-m)- 
In the limit L —> oo, the ratio fj^9(L\&;r,p)/L approaches the one-loop term of 
the reduced bulk free-energy density /b(<j;f, p), namely 

n f (d) 

/b,[i](o-;T,p) = - / ]nn & (q\r,p). (C.3) 
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Subtracting from equation (jC.ip the bulk term and using Poisson's summation 
formula (|5.58[> gives 

/r™± ;[ i] r, P) = fZffi(L\&; f, p) - Lf h>[1] (a; f, p) 

oo Ad) 

= n ^2 cos(q m jL)\nfla{q\T,p). (C.4) 
j=i U 

The summands in the second line of this equation may be conveniently expressed 
in terms of the free bulk propagator G^ +2,m \ To show this, we insert the identity 
1 = V p • p/(d — m) into the integral J^ d and integrate by parts, obtaining 
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where P is a (d + 2 — m)-dimensional wave vector. We have exploited the facts that 
Oo-(g|r, p) depends on p merely via p 2 and that p-integrals over rotationally invariant 
functions h(p 2 ) simplify to 

AD) poo 

/ h{p 2 ) = K D / dpp^hip 2 ), (C.6) 
Jp Jo 

with 

K D = 2 1 - d tt- d ' 2 /T{D/2), (C.7) 

the standard geometrical factor resulting from the angular integrations. 
The result (|C5|) suggests to introduce the functions 

r(d) 2cos(g m j) 



P r ^ = f 2cos( gm j) 
mA,P) -f^J q T + pk 2 +p 2 + ki 

oo 

= 2^G^(^|l;r,p) ) (C.8) 



j'=i 



where G^' m \x\a; f , p) is the free bulk propagator (|4. 1|) . It is straightforward to show 
that these functions have the property 

^P m .d+2(.T, P ) =-±-P m , d ( T ,p), (C.9) 

which they share with the functions Qd,i{v)/y used in [3^ and [IU] [cf. equations (4.17) 
and (B5) of the first and second reference, respectively]. Let us differentiate the series 
of P m ,d+2 termwise with respect to r and interchange the differentiation with the 

momentum integration ■ The r-derivative of the integrand can be replaced by its 
derivative with respect to p 2 . An integration by parts as in equation (IC5|) then gives 
relation (lCJ]) . 

Upon combining the above results, we finally arrive at 

/™±;[i]M*;^) = ^Trna^^ 2 L 1 - 2 ^ P m , d+2 (fL 4 /a,pL 2 /a). (CIO) 

In order to determine the Casimir amplitudes A^ BC (d, m, n), we need the value 
of P m ,d+2 at f = p = 0, which we denote as P^ d+2 - To compute P^ d , we substitute 
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the scaling form (I4.3[) of the free bulk propagator into equation (|C.8I) and use the 
asymptotic behavior (I4.8[) of the scaling function & my d, obtaining 

oo 

PZ = *W(0,0) = 2j2Gl d ' m) (jz\l;0,0) 
i=i 

= 2$ I ( ^C(2d-m-4), (C.ll) 

where is the coefhcient defined in equation (14.91) . Its insertion into 

equation (jC.lip yields the expression quoted in equation (I5.47[) . From equation (|C.10|) 
one can read off the expression for the Gaussian Casimir amplitude A^q given in 
the first line of equation (|5.46[) . The result presented in its second line follows upon 
insertion of the result for P^ d+2 implied by equation (|5.47p . 
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